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ABSTRACT 

on : In this first paper of a series, we report on the use of quasar spectra obtained with the 

UVES high resolution spectrograph and available through the ESO VLT archive to 
build the first homogeneous and unbiased sample of sub-DLA systems, absorbers with 
H I column densities > 10 19 cm ~ 2 but lower than the classical definition of damped 
Lya systems (DLAs) 2 x 10 20 cm" 2 . A systematic investigation of the properties 
, of these systems and a comparison with those of the DLAs (Paper II of this series; 

Peroux et al. 2003b) is expected to provide new clues on the association of high column 
density absorbers with galaxies and on the overall evolution of the neutral hydrogen 
gas mass and metal content in the Universe. In the spectra of 22 quasars which were 
found suitable for a sub-DLA search, we identified 12 sub-DLAs and 1 borderline 
Jjj | case between the DLA and sub-DLA systems in the redshift interval z = 1.8 — 4.3. 

We measured the column densities of Hi and of up to 16 additional ions of low-, 
intermediate- and high-ionization, Oi, Cn, Civ, Sin, Siiv, Ni, Sn, Mgl, Mgn, Alii, 
Aim, Fen, Fein, Nin, Znii, and Cm. We further investigated the significance of 
the ionization corrections in the determination of the chemical abundances from the 
low-ionization ions in the sub-DLA H I column density range. Using the predictions 
of different ion ratios as a function of the ionization parameter computed with the 
CLOUDY software package, we have estimated that with the exception of one case, 
the ionization corrections to the abundances of 9 systems for which we were able 
to constrain the ionization parameter, are lower than 0.2 dex for all of the elements 
except Alii and Znii down to Hi column densities of logA(Hi) = 19.3 cm" 2 . We 
finally present the first sub-DLA chemical abundance database which contains the 
abundance measurements of 11 different elements (O, C, Si, N, S, Mg, Al, Fe, Ni, Zn, 
and Cr). We took advantage of the lower Hi column densities in sub-DLAs to measure, 
in particular, the O and C abundances using lines which are normally saturated in 
DLAs. 

Key words: galaxies: abundances - galaxies: high-redshift - quasars: absorption 
lines - quasars. 
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1 INTRODUCTION 

Galaxy formation and evolution remain important issues in 
our understanding of the early epoch of the Universe. Funda- 
mental steps towards these issues would be the reconstruc- 
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tion of the chemical histories of galaxies and the metal cen- 
sus of the Universe. They are investigated both via the tra- 
ditional approach of observing the galaxy starlight, and via 
the absorption line systems detected along the quasar lines 
of sight independently of their distance, luminosity and mor- 
phology. The absorption line systems, in particular, provide 
a very accurate observational method to measure the neu- 
tral gas and metal content of the Universe up to very high 
redshifts, z > 4 (e.g. Lu et al. 1996a; Prochaska, Gawiser 
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& Wolfe 2001; Dessauges-Zavadsky et al. 2001a; Songaila & 
Cowie 2002). 

The quasar absorption line systems are divided into 
three classes according to their neutral hydrogen column 
density: the Lya forest with iV(Hl) ranging from ~ 10 12 to 
1.6 x 10 17 cm -2 , the Lyman limit systems (LLs) with iV(H i) 
> 1.6 x 10 17 cm -2 , and the damped Lya systems (DLAs) 
with JV(Hl) > 2.0 x 10 20 cm" 2 (Wolfe et al. 1986). They 
thus probe different media from voids to halos and disks of 
both dwarf and normal (proto-)galaxies (Prochaska & Wolfe 
1998; Ledoux et al. 1998; Haehnelt, Steinmetz & Rauch 
1998; Matteucci, Molaro & Vladilo 1997; Jimenez, Bowen 
& Matteucci 1999; Calura, Matteucci & Vladilo 2003). The 
DLAs, in particular, are a cosmologically important pop- 
ulation, since they dominate the neutral hydrogen content 
of the Universe (e.g. Wolfe et al. 1995; Storrie-Lombardi & 
Wolfe 2000) and are believed to be the progenitors of the 
present-day galaxies. 

A recent study by Peroux et al. (2003a) has, however, 
suggested that the relative importance of lower H I column 
density systems, the so-called 'sub-DLA systems' with H I 
column densities between 10 19 and 2 x 10 20 cm" 2 , increases 
with redshift. Indeed, this until now poorly studied class of 
absorbers should contain about 45% of the neutral hydrogen 
mass at z > 3.5, and may thus play an important role at 
high redshift in the study of the redshift evolutions of both 
the total amount of neutral gas and the metal content in the 
Universe. The sub-DLAs, similarly to the DLAs, present the 
advantage that the H I column density of these systems can 
accurately be estimated thanks to the presence of damp- 
ing wings in their Lya lines, observed down to Hi column 
densities of 10 19 cm -2 . 

At the present time only the chemical abundances of 
DLAs have been used to trace the metal content in the Uni- 
verse, but this has shown unexpected results. Contrary to 
the cosmic chemical evolution model predictions (e.g Pei & 
Fall 1995; Pei, Fall & Hauser 1999; Cen & Ostriker 1999; 
Cen et al. 2003), the most recent observations indicate only 
a mild evolution of the metal content with redshift (Pettini 
et al. 1999; Prochaska & Wolfe 2000; Vladilo et al. 2000; 
Dessauges-Zavadsky et al. 2001b; Prochaska & Wolfe 2002). 
Only Savaglio (2000), who focused on lower Hi column den- 
sity DLA systems, claims evidence for redshift evolution. 
If the Peroux et al. (2003a) predictions are confirmed, the 
sub-DLA metallicities should be taken into account in or- 
der to obtain a complete picture of the redshift evolution of 
the Universe's metal content. Indeed, the current metallicity 
studies which focused only on the higher column density ab- 
sorption line systems, may have provided a biased and/or in- 
complete view of the global metal evolution in the Universe 
at z > 3.5. The observational confirmation of the impact 
of the sub-DLA systems at high redshift thus is crucial. In 
addition, the investigations of the chemical, kinematic and 
clustering properties of the sub-DLAs to see whether we are 
dealing with a different class of objects with a lower H I 
content /mass relative to the DLAs are our further motiva- 
tions for studying these absorbers. If sub-DLAs are the basic 
building blocks of the hierarchical growth of galactic struc- 
tures as suggested by Peroux et al. (2003a), their abundance 
analysis may provide an insight into the very early stages of 
the chemical enrichment of galaxies. 

The work presented in this paper therefore aims at con- 



structing the first homogeneous sample of sub-DLA systems. 
For this purpose, we have used all of the high-resolution 
quasar spectra available to us in July 2001 in the ESO 
UVES-VLT archives. Section 2 explains how the quasar 
spectra were selected and how the sub-DLA systems were 
identified. We provide a detailed chemical analysis for each 
sub-DLA system in Section 3, and we address the issues of 
the photoionization effects in Section 4. Finally, in Section 5, 
we summarize all the ionic column density and the absolute 
abundance measurements. In Paper II of this series (Peroux 
et al. 2003b), we use the database presented here in combi- 
nation with the data on DLAs and other sub-DLAs from the 
literature to discuss the overall properties of these objects. 



2 CONSTRUCTION OF A HOMOGENEOUS 
SUB-DLA SAMPLE 

All the observations presented in this paper were acquired 
with the Ultraviolet- Visual Echelle Spectrograph (UVES; 
see D'Odorico et al. 2000) on the VLT 8.2 m ESO telescope 
at Paranal, Chile, between the Commissioning and the Sci- 
ence Verification time of the instrument at the end of 1999 
and beginning of 2000, respectively, and July 2000. All these 
data reached in July 2001 the end of their one-year long pro- 
priety period and could be retrieved from the ESO archives. 
A total of 35 high resolution, high signal-to-noise ratio (S /N) 
quasar spectra was available, illustrating the success of both 
the instrument and the data distribution policy of the ESO 
archives. 

2.1 Quasar Selection 

The goal of our study being well defined, we have selected 
from the 35 quasars available in the ESO UVES archives 
only the quasars suitable for the sub-DLA search. First, we 
took the conservative approach of excluding all quasars ex- 
hibiting Broad Absorption Lines (BALs), since the nature of 
these objects and their characteristic features are not fully 
understood yet, rendering the distinction between true in- 
tervening absorption systems and other quasar spectral fea- 
tures difficult. However, we have not excluded the quasars 
selected for the study of associated systems, z a b s — z cra , 
since the quasars considered in the ESO programme 65. P- 
0038 (see Table do not exhibit BALs which are often 
confused with associated absorbers. Secondly, we have ex- 
cluded the quasars which were observed for the sole purpose 
of studying the Lya forest (e.g. for the analysis of voids, 
low column density distribution) . Indeed, these quasars have 
been pre-selected for having no high column density systems 
in their spectrum, precisely the type of systems we are look- 
ing for in the present study. In one case, however, a sub-DLA 
has been found in the later sample of quasars, the quasar 
Q1101-264 at z cm = 2.12 observed during the UVES Sci- 
ence Verification. This absorber is included in the present 
abundance analysis, but for consistency reasons it will be 
excluded from the study of the sub-DLA statistical proper- 
ties (see Paper II). 

The above selection process resulted in a sample of 22 
quasars. A large majority of these quasars were observed 
with UVES because of their intervening DLA systems. There 
is a priori no reason why this should create a bias in our 
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Table 1. Observational properties of the selected quasar spectra from the ESO UVES archives and their intervening DLAs and sub-DLAs 



Quasar name z 6m ESO prog. No. Settings Wavelength coverage* Exposure zdla 2 sub— DLA 

A s +(Ref.) 



Q0000-2620 


4.11 


Commissioning 


B437+R860 


3770- 


-5000, 6700- 


-10500 


19700 


3.390 




APM BR J0307— 4945 


4.75 


Commissioning 


R600, R760, R800 


5000- 


-9900 




26300 


4.466 






3.23 


Commissioning 


R4274-R860 


3655- 


-4888 6722- 


-10000 


9500 


3.025 




>^uu^i J. ill? 


2.50 


65 O-0063 


R41 2+R860 


3550- 


-4740 6710- 

t: 1 5 U 1 J.VJ 


-10000 


5400 


2.375 
2.476 




HE 094(1—1050 


3.05 


65 O-0474 


R34fi+R "SSO 


3600- 


-3870 4780- 


-6809 


3600 






Q1038-272 


2.32 


65.O-0063 


B437+R860 


3765- 


-4998, 6707- 


-10000 


5400 






Q1101-264 b 


2.14 


Science Verification 


B346+R580, B437+R860 


3050- 


-8530 




45000 




1.838 6 (1) 




2.31 




R^fiD Rfi70 
douU) nui u 


3184- 


— £ ±uzr i ± ) ouyo 


-7714 


10880 


1.662 




Q1151+068 


2.76 


65.O-0158 


B346+R580, B437+R860 


3100- 


-10000 




21600 


1.774 




PKS 1 1 57+01 4 


1.99 


65 O-0063 


R34fi+R ^0 


3050- 


-3870 4500- 


-6500 


7200 


1.944 




Q1223+1753 


2.94 


65.P-0038, 65.O-0158 


R390+R564, R437+R860 


3300- 


-10000 




18000 


2.466 


2.557 


PKS 1232+0815 


2.57 


65.P-0038 


B390+R564 


3300- 


-6650 




10800 


2.338 




Q1409+095 


2.86 


65.O-0158 


B437+R860 


3769- 


-4998, 6707- 


-10000 


18000 


2.456 


2.668 (2,3,4) 


Q1444+014 


2.21 


65.O-0158 


B390+R564 


3294- 


-6650 




18000 




2.087 (3,5) 


Q1451+123 


3.25 


65.O-0063 


B424+R750 


3640- 


-4855, 5900- 


-9400 


7200 


2.469 


2.255 d (2) 
3.171 e (1,2) 


Q1511+090 


2.88 


65.P-0038 


B390+R564, B437+R860 


3300- 


-10000 




32400 




2.088 


Q2059-360 


3.09 


65.O-0063 


B480+R740 


4200- 


-5340, 5900- 


-9300 


10800 


3.083 


2.507 


Q2116-358 


2.34 


65.O-0158 


B390+R564 


3294- 


-6650 




7200 




1.996 (6,7,8) 


Q2132-433 


2.42 


65.O-0158 


B390+R564 


3294- 


-6650 




3600 


1.914 




PSS J2155+1358 


4.26 


65.O-0296 


R565, R800 


4633- 


-9800 




11600 


3.316 


3.142 
3.565 
4.212 


Q2206-199 


2.56 


65.O-0158 


B390+R564, B437+R860 


3294- 


-10000 




25200 


1.920 
2.076 




PSS J2344+0342 


4.24 


65.O-0296 


R565, R800 


4633- 


-9800 




13500 


3.220 


3.882 



* With some gaps in the wavelength coverage. 

a The detection of sub-DLAs is uncertain (see Section l3.ll . 

b This quasar has been observed for the study of the Lyct forest. Nevertheless, an analysis of the spectra 
revealed the presence of a sub-DLA. For consistency reasons, this system will not be included in the study 
of the sub-DLA statistical properties presented in Paper II. 

c This quasar is gravitationally lensed. Only the brightest line of sight was included in our study. 

d This absorption line system is a borderline case between the DLA and sub-DLA systems with a H I column 

density of log iV(Hl) = 20.30 ± 0.15. It is not included in the sub-DLA sample. 

e The spectra do not cover the Lya absorption line at A t, s ~ 5070 A of this sub-DLA. It has been identified 
thanks to the Hi column density measurement reported by Petitjean, Srianand Ik, Ledoux (2000). 
References to the sub-DLAs already reported in the literature: 

(1) Petitjean, Srianand & Ledoux (2000); (2) Lanzetta et al. (1991); (3) Pettini et al. (2002); (4) Lu et al. 
(1993); (5) Wolfe et al. (1995); (6) Lanzetta, Wolfe k, Turnshek (1987); (7) Wampler, Bergeron & Petitjean 
(1993); (8) M0ller, Jakobsen & Perryman (1994) 
Observing programmes: 

65.P-0038: "Probing the nuclear regions of QSOs with z^bs — z cm associated systems" 

(Srianand, Petitjean & Aracil) 
65.O-0063: "Molecular hydrogen at high redshift" (Ledoux, Petitjean, Srianand & Rauch) 
65.O-0158: "Abundance ratios in Damped Lya systems: clues to the nucleosynthesis of N and Mn" 

(Pettini, Bergeron & Petitjean) 
65.O-0296: "A study of the IGM-galaxy connection at z = 3 using the UVES high quality absorption spectra" 

(D'Odorico, Cristiani, Fontana, Giallongo & Savaglio) 
65.O-0474: "Deuterium abundance in the high redshift QSO absorption system towards HE 0940—1050" 

(Molaro, Centurion & Bonifacio) 



sample. We assume that the sub-DLA systems are not pref- 
erentially found close to the DLA systems. An analysis of 
the clustering properties of the sub-DLA systems identified 
in our sample is presented in Paper II. 

Table Q summarizes the observing programme details 
(i.e. the ESO programme number, the exposure time, the 
settings), the wavelength coverages, the redshifts of the 22 
selected quasars and the redshifts of their intervening DLA 
systems. The sample of the selected quasars should result in 



the construction of an unbiased and homogeneous sample of 
sub-DLA systems. 

2.2 Data Reduction 

The constructed sample of 22 quasars has benefited from 
a homogeneous data reduction. All of the quasar spectra 
were reduced using the ESO data reduction package MIDAS 
and the UVES pipeline in an interactive mode available as a 
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Figure 1. Output of the automated sub-DLA search algorithm for the quasar Q2116— 358. The solid line corresponds to the equivalent 
width spectrum and the dotted line to the error array. The straight dashed line represents the threshold above which the sub-DLA 
systems are detected. In this case a sub-DLA at ~ 3640 A (z a b s = 1.996) is clearly identified. 



MIDAS context. A detailed description of the pipeline can be 
found in Ballester et al. (2000). To optimize the results we 
made a systematic check of each step of the pipeline reduc- 
tion. Once reduced, the wavelengths of the one-dimensional 
spectra were converted to a vacuum-heliocentric scale, the 
individual spectra were co-added using their S/N ratio as 
weights to get the maximum signal-to-noise ratio, and the 
various settings were combined to obtain the maximum 
wavelength spectral coverage. An average resolution of 5 — 7 
km s~ x FWHM was achieved in the spectra, and an average 
S/N per pixel between 10 and 150 is observed from spec- 
trum to spectrum according to the exposure time and the 
magnitude of the quasar. 

The final step was the normalization of the resulting 
spectra obtained by dividing them by a spline function fit- 
ted to smoothly connect the regions free from absorption 
features. The continuum in the Lya forest was fitted by 
using small regions deemed to be free of absorptions and 
by interpolating between these regions with a spline. The 
higher the redshift of the quasar, the more difficult is the 
normalization in the Lya forest due to the higher density of 
absorption lines (e.g. Kim, Cristiani & D'Odorico 2001). The 
procedure we applied for the normalization is the same as 
the one we used for two high-redshift and four intermediate- 
redshift quasars we have previously analysed for a detailed 
chemical abundance study of their intervening DLA systems 
(Dessauges-Zavadsky et al. 2001a; Levshakov et al. 2002; 
Dessauges-Zavadsky et al. 2002). This procedure has proven 
to lead to an accuracy of the normalization level of about 
5 — 10% in the Lya forest and better than 5% redwards of the 
Lya emission of the quasar (according to the signal-to-noise 
ratio of the spectra). For each quasar we tested the derived 
normalization first by selecting different regions in the con- 
tinuum and by varying the number of these regions used for 
the fit of the continuum with a spline. In the Lya forest, in 
particular, some regions are key regions in the continuum 
definition, when omitted, the shape and level of the con- 



tinuum are drastically modified. We tried to qualitatively 
and quantitatively estimate these modifications. Secondly, 
we checked whether the adopted continuum fit is within the 
1 a error of the spectra along the whole wavelength range 
sampled, and we estimated the normalization accuracy by 
determining the mean deviation of the adopted continuum 
fit from the observed continuum within its 1 a error and the 
maximum variation of this fit tolerated within the 1 a error 
of the spectra. 

2.3 Sub-DLA Identification 

The sample of quasars to be analysed being defined and their 
high resolution spectra being reduced, we used these spec- 
tra to identify their intervening sub-DLA systems. In order 
to select all the sub-DLAs, we have used a detection algo- 
rithm equivalent to the one applied by Peroux et al. (2001) 
to detect the DLAs in medium resolution spectra (~ 2 A 
pixel -1 ), and which was previously used by Lanzetta et al. 
(1991), Wolfe et al. (1995), and Storrie-Lombardi & Wolfe 
(2000). This technique consists in building an equivalent 
width spectrum over 500 pixel wide boxes for each quasar. 
The analysis was done 3000 km s~ x bluewards from the Lya 
emission of the quasar in order to avoid possible contami- 
nations from systems associated with the quasar itself, and 
it ended where the S/N was too low to detect absorption 
features at the sub-DLA threshold, i.e. EW rcs t = 2.5 A. 

Fig.^displays an example of the output of the sub-DLA 
search algorithm in the case of the quasar Q2116— 358. The 
solid line is the equivalent width spectrum and the dotted 
line is the error array. The straight dashed line corresponds 
to the threshold above which the sub-DLA systems are de- 
tected. In this example, the sub-DLA system at ~ 3640 A 
(.Jabs = 1.996) is clearly identified. 

We then supplemented the sub-DLA candidates found 
by the algorithm with a visual identification. Indeed, since 
the algorithm selects any region with the required equivalent 
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Figure 2. Lyct absorption line profiles for four different Hi col- 
umn densities in the sub-DLA range: from left to right, log ./V(H i) 
= 20.3 (solid line), 20.0 (dotted line), 19.5 (solid line) and 19.0 
(dotted line) with the Doppler parameter b fixed at 20 km s _1 . 
This illustrates that in the sub-DLA range, 10 19 < AT(Hl) 
< 2 X 10 20 cm -2 , the damping wings are clearly observed. 

width regardless of the shape of the absorption feature and 
whether or not the absorption goes to zero, an additional 
check of the candidates is necessary to distinguish the true 
sub-DLAs from blends of several Lya clouds. Three main vi- 
sual criteria were used to identify a sub-DLA system: i) the 
saturation of the absorption feature identified as the Lya ab- 
sorption line of the sub-DLA candidate, ii) the presence of 
damping wings in the Lya absorption lines which are clearly 
observed in the sub-DLA H I column density regime between 
10 19 and 2 x 10 20 cm" 2 as illustrated in Fig.H and Hi) the 
detection of lines of higher members of the Lyman series 
with the same H I column density as the one measured from 
the Lya line when the wavelength coverage is available. An 
additional visual criterion, though not necessary, was the 
identification of associated metal lines which are very likely 
to occur in absorbers with H I column densities greater than 
10 19 cm" 2 . On the basis of these visual criteria and to im- 
partially select the sub-DLAs, three of us (M.D.-Z., CP., 
T.-S.K.) looked independently at the sub-DLA candidates 
found by the algorithm and closely inspected the original 
spectra to determine whether or not these absorption fea- 
tures correspond to sub-DLAs. 

We reached a good agreement in all but one case, the 
quasar APM BR J0307-4945 at z em = 4.78. Given the high 
redshift of this quasar the Lya forest is very dense, and 
thus at three different redshifts the presence of sub-DLA 
systems is ambiguous (see Section r3.ll . The absence of lines 
of higher members of the Lyman series for two of them does 
not allow us to draw a conclusion on the true Hi column 
density of these systems, leading us to not include them in 
the sub-DLA sample. In the third candidate, the Lyman 
lines are heavily blended and prevent us from obtaining an 
accurate Hi column density measurement, thus this candi- 
date is also not analysed and not included in our sub-DLA 
sample. Moreover, the sub-DLA system at z a b s = 3.171 to- 
wards Q1451+123 has been identified thanks to the Hi col- 
umn density measurement reported by Petitjean, Srianand 
& Ledoux (2000), the spectra we had at our disposal do 
not cover the Lya line of this system. With the lines of 
higher members of the Lyman series — Ly/3, Ly7 and Ly8 — 
we confirm the reported H I column density measurement of 
this absorption system (see Section r3.6.21 . This system is in- 
cluded in the abundance analysis, but it is excluded from the 
study of the sub-DLA statistical properties (see Paper II), 
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Figure 3. Hi column density distribution of the 12 sub-DLA 
systems and 1 borderline case between the DLA and sub-DLA 
systems detected towards the 22 selected quasars. 

similarly to the sub-DLA system towards Q1101 — 264 (see 
the explanation in Section 12.111 . Finally, the identified ab- 
sorption line system at z a b s = 2.255 towards Q1451+123 is 
a borderline case between the DLA and sub-DLA systems 
with a Hi column density of log A r (Hl) = 20.3 ± 0.15 (see 
Section f3.6.H . We provide a detailed abundance analysis of 
this system, but we do not consider it as a sub-DLA system, 
and thus we do not include it in the sub-DLA sample used 
for the different studies undertaken in Paper II. 

In Table^we give the redshifts and the references of the 
12 sub-DLAs and 1 borderline case between the DLA and 
sub-DLA systems detected towards the 22 selected quasars, 
and in Fig. [3]we show their H I column density distribution. 
The detected systems span a range of Hi column densities 
from 10 19 ' 3 to 10 20 3 cm . No particular trend can be con- 
firmed with this small number of objects yet. The sample 
of the 10 sub-DLAs, when excluding the sub-DLA towards 
Q1101— 264 and towards Q1451+123, constitutes an unbi- 
ased and homogeneous sample of sub-DLAs, ideal for a sta- 
tistical analysis (see Paper II). 



3 IONIC COLUMN DENSITIES 

In this Section we present the ionic column density mea- 
surements of the 12 detected sub-DLA systems and 1 bor- 
derline case between the DLA and sub-DLA systems. All of 
the column densities of the metal species were derived with 
the Voigt profile fitting technique. This technique consists in 
fitting theoretical Voigt profiles to the observed absorption 
metal lines well described as a complex of components, each 
defined by a redshift z, a Doppler parameter b, a column 
density and the corresponding errors. The fits were per- 
formed using a x 2 minimization routine f itlyman in MIDAS 
(Fontana & Ballester 1995). 

We assumed that all of the metal species with similar 
ionization potentials can be fitted using identical component 
fitting parameters, i.e. the same b (which means that macro- 
turbulent motions dominate over thermal broadening) and 
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Table 2. Component structure of the z a t, B = 1-838 sub-DLA system towards Q1101— 264 



No 2 abs v* el 


b(a b ) 


Ion 


logW((T logJV ) 


No 


■^abs 


^rcl 


b(a b ) 


Ion 


logiV(<Tl og jv) 


km s 


km s — 1 




cm -2 






km s _1 


km s _1 




_9 

cm 


Low-ion transitions 


1 1.837407 -158 


3.0(0.9) 


Cn 


12.30(0.03) 


8 


1.838902 





5.3(0.1) 


Sin 
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the same z in the same component, and allowing for vari- 
ations from metal species to metal species in N only. We 
distinguish three categories of metal species with similar 
ionization potentials: the low-ion transitions (i.e. the neu- 
tral and singly ionized species), the intermediate-ion tran- 
sitions (e.g. Fein, Aim) and the high-ion transitions (e.g. 
Civ, Siiv). By using relatively strong lines (but not satu- 
rated) to fix the component fitting parameters (the b and z 
values for each component), we then obtain excellent fit re- 



sults even for weak metal lines and for metal lines in the Lya 
forest where the probability of blending is high, by allowing 
only the column density to vary. For all of the analysed sys- 
tems exhibiting multicomponent velocity structures, we had 
a sufficient number of relatively strong metal lines to con- 
strain well the fitting parameters. 

Tables 01— 1 1 41 present the results of the component per 
component ionic column density measurements for the fit- 
ting model solutions of the low-, intermediate- and high-ion 
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transitions for all of the analysed systems. The reported er- 
rors are the 1 a errors computed by f itlyman. These errors 
do not take into account the uncertainties on the continuum 
level determination. For the saturated components, the col- 
umn densities are listed as lower limits. The values reported 
as upper limits correspond to cases contaminated by signif- 
icant line blendings due to Hi clouds in the Lya forest or 
telluric lines. Figs. |5T- |2"§I (odd numbers) show the best fitting 
solutions of the low-, intermediate- and high-ion transitions 
for all of the analysed systems. In these velocity plots, v — 
corresponds to an arbitrary component, and all the identified 
components are marked by small vertical bars. The thin solid 
line represents the best fit solution. We point out that the 
high-ion transitions in the bulk of sub-DLA systems present 
very different profiles from the low-ion transitions, whereas 
the intermediate-ion transitions show relatively similar pro- 
files with small differences in the model fitting parameters. 
This will be discussed individually for each system in the 
following sub-sections. 

The neutral hydrogen column densities were measured 
from the fits of the Lya damping line profiles. The fe-value 
was usually fixed at 20 km s~ or left as a free parameter 
— similarly to the DLA systems the fit of the damped Lya 
line of sub-DLA systems is independent of the b- value — , 
and the redshift z was fixed at the redshift of one of the 
strongest components of the low-ion metal line profiles or 
also left as a free parameter depending on the sub-DLA sys- 
tem (see the comments in the following sub-sections). When 
other lines of higher members of the Lyman series were cov- 
ered in the spectra, we used them to check the H I column 
density obtained from the fit of the Lya line. Figs. IH- I28I 
(even numbers) show the results of the Hi fitting solutions 
for all of the analysed systems. The small vertical bar cor- 
responds to the redshift used in the best fit solution and 
the thin solid line represents the best fit solution. In one 
case, the sub-DLA at z a b s = 3.171 towards Q1451+123, the 
Lya line is not covered in the available spectra, thus the H I 
column density measurement was performed using the Ly/3, 
Ly7 and Ly8 lines, the other members of the Lyman series 
being heavily blended. 

Throughout the analysis we have adopted the list of 
the atomic data — laboratory wavelengths and oscillator 
strengths — compiled by J. X. Prochaska and collaborators, 
and presented on the web site "The HIRES Damped Lya 
Abundance Database . The most recent measurements of 
the A and /-values of the metal-ions that impact the abun- 
dances of high-redshift absorption systems and their refer- 
ences are reported there. 

We now comment on the individual sub-DLA systems. 

3.1 APM BR J0307 4945 

The spectra of this quasar (Storrie-Lombardi et al. 2001) 
as well as the chemical abundances of one of the high- 
est redshifted damped Lya system discovered by Peroux 
et al. (2001) have been analyzed in detail by Dessauges- 
Zavadsky et al. (2001a). Given the high redshift of this 
quasar (z om = 4.75) several intervening sub-DLA systems 
are expected. Three candidates at z^s = 3.356, 3.589 and 

1 http:/ /kingpin. ucsd.edu/~hiresdla/ 
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Figure 4. Normalized UVES spectrum of Q1101— 264 showing 
the sub-DLA Lya profile with the Voigt profile fit. The zero ve- 
locity is fixed at z = 1.838902. The vertical bar corresponds to 
the velocity centroid used for the best fit, z = 1.838608. The 
measured H I column density is log ./V(H i) = 19.50 ± 0.05. 



4.214 have been detected. However, the high density of ab- 
sorption lines in the Lya forest of such distant objects makes 
the blends of several absorption clouds very likely to occur. 
The lines of higher members of the Lyman series are nec- 
essary to distinguish blends of clouds from high Hi column 
density absorbers. As discussed below, none of the three sub- 
DLA candidates can be confirmed in the considered spectra 
and they probably never will be. They thus are not included 
in any of our further analysis (this paper and Paper II). 

For the two candidates at lower redshifts, z a hs ~ 3.356 
and 3.589, we do not have any information on the higher 
members of the Lyman series, these lines being beyond the 
quasar flux cut-off. Moreover, for the system at z a b s = 3.356 
no associated metal lines are observed. Thus, since the Lya 
absorption lines of these two candidates do not show clearly 
any damping wing and since they can be fitted with one 
single system with iV(Hl) > 10 19 cm" 2 as well as with sev- 
eral clouds of lower Hi column densities, the identification 
of these sub-DLA candidates is doubtful, and we choose to 
consider these two detected absorption line systems as un- 
likely sub-DLAs. 

Several metal lines (Cn, Cm, Civ, Sin, Sim, Siiv, 
Nm, Aim) and the hydrogen lines Lya, Ly/3 and Ly7 have 
been identified in the candidate at z a b s = 4.214. From the 
metal line profiles two main sub-systems can be clearly dis- 
tinguished, separated by ~ 400 km s _1 . These two sub- 
systems have already been carefully analysed by Levshakov 
et al. (2003). From their Monte Carlo Inversion analysis 
they determined the Hi column densities of 1.4 x 10 17 and 
2.3 x 10 17 cm -2 in the bluest and reddest sub-systems, re- 
spectively, corresponding to values characterizing the class 
of the Lyman limit absorbers. Although these recovered val- 
ues are self-consistent, they are not necessarily a unique so- 
lution. Indeed, the iV(Hl) measurement in the reddest sys- 
tem is well constrained by the three observed hydrogen lines. 
However, all of the three hydrogen lines in the bluest sys- 
tem are highly blended, hence the real AT(Hl) value may 
be higher and can possibly be above the sub-DLA AT (Hi) 
definition threshold. Given this uncertainty, we prefer again 
to consider this absorption line system only as a sub-DLA 
candidate. 
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Figure 5. Velocity plots of the low-, intermediate- and high-ion transitions (normalized intensities shown by dots with 1 a error bars) for 
the sub-DLA towards Q1101— 264. The zero velocity is fixed at z = 1.838902. The vertical bars mark the positions of the fitted velocity 
components (see Table |2"1. For this and all the following figures, the notation "(Lya)" near the name of a transition means that the line 
is located in the Lya forest, and the symbols © correspond to telluric lines. 
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Table 3. Component structure of the z a t, B = 2.557 sub-DLA system towards Q1223+1753 
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3.2 Q1101-264, a sub _DLA = 1-838 

This quasar has been observed as part of the UVES Science 
Verification programme for the study of the Lya forest (see 
Kim, Cristiani & D'Odorico 2002). It has previously been 
observed and analysed at lower resolution, and the absorp- 
tion line system at 2 a b s = 1.838 we are interested in has 
already been identified (e.g. Carswell et al. 1984). Thanks 
to the large amount of UVES observing time devoted to this 
quasar, the spectra show an extremely high S/N. The qual- 
ity of the data allowed us to measure with high accuracy 
the column densities of 11 ions and 7 elements (Si, O, S, 
C, Al, Fe, Mg) in the sub-DLA system at z abs = 1.838 ob- 
served in this line-of-sight. 4 a upper limits on the column 
densities of Znn and Cm have been provided, since these 



lines are not detected in the spectra. Only the hydrogen line, 
Lya, is observed in the available spectral coverage, but the 
absence of contamination and the presence of well defined 
damping wings lead to a very reliable JV(H i) measurement 
of 19.50 ± 0.05 obtained by leaving the redshift and the b- 
value as free parameters (see Fig. 2J. This absorber has al- 
ready been studied by Petitjean, Srianand & Ledoux (2000) 
with regards to its molecular content. We confirm their Hi 
column density measurement and their [Fe/H] 2 and [S/H] 
abundance measurements. 

The component structure of this absorber is presented 



2 [X/H] = log[Ar(X)/jV(H)] DLA -log[7V(X)/AT(H)] . 
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Figure 7. Same as Fig.0for the sub-DLA towards Q1223+1753. The zero velocity is fixed at z = 2.557640. The vertical bars mark the 
positions of the fitted velocity components (see Table 151 . 
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Figure 6. Normalized UVES spectrum of Q1223+1753 showing 
the sub-DLA Lya profile with the Voigt profile fit. The zero veloc- 
ity is fixed at z = 2.557640. The vertical bar corresponds to the 
velocity centroid of the strongest component of the low-ion tran- 
sition lines used for the best fit, the component 6 at z = 2.557640. 
The measured H I column density is log N(H i) = 19.32 ± 0.15. 



km s~ , while the high-ion transitions are well fitted with 
only 6 components (see Fig. |SJ . This difference in the line 
profiles suggests that the low- and high-ion lines are not 
emitted from the same physical region, which is a commonly 
observed situation in DLAs. Nevertheless, in this partic- 
ular case, surprisingly the lines of C +++ and Si +++ , two 
ions with the same ionization degree, do not show the same 
line profiles and require slightly different fitting parameters 
(z, b and N). In addition, the Aim lines better resemble 
the high-ion transition profiles than the low-ion ones, con- 
trary to what is usually observed in DLA studies where the 
intermediate-ion line profiles are very similar to the low- 
ion line profiles, indicating they are emitted from the same 
physical region. Unfortunately, no other intermediate-ion is 
detected in this sub-DLA system to confirm these observa- 
tions. 



in Table |21 for the low-, intermediate- and high-ion transi- 
tions. The low-ion profiles are described with a relatively 
complex structure of 11 components extended over ~ 210 



3.3 Q1223+1753, z sub - D LA = 2.557 

This quasar, also called LBQS 1223+1753, was discovered 
by Foltz et al. (1987) and it has been observed on many 



A Homogeneous Sample of Sub-DLAs I: Sample and Chemical Abundance Measurements 11 



Table 4. Component structure of the z abs = 2.668 sub-DLA 
system towards Q1409+095 
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13.32(0.04) 




Oi 


14.24(0.15) 




Cn 


> 16.00 




Fen 


13.63(0.24) 


2 2.668227 4.8(0.4) 


Si II 


14.06(0.03) 




Oi 


15.27(0.11) 




Cn 


> 15.37 




Sn 


13.54(0.06) 




Ni 


< 13.49 




Cm 


12.49(0.15) 




Fen 


13.79(0.05) 




Fe in 


13.17(0.12) 




Aim 


11.89(0.05) 



* Velocity relative to z = 2.668227 



occasions for the study of its intervening DLA system at 
Zabs = 2.466. The large H I column density of this absorber 
allows the analysis of a great number of transitions (over 
20) as well as of the H2 molecule (e.g. Pettini et al. 1994; 
Centurion et al. 2000; Petitjean, Srianand & Ledoux 2000; 
Ellison, Ryan & Prochaska 2001; Prochaska et al. 2001). 
In addition to this DLA, we identified a sub-DLA system 
at z abs = 2.557 with a logTV(Hl) of 19.32 ± 0.15, and we 
present here for the first time its chemical analysis. 

This sub-DLA system was identified thanks to its Lya 
absorption line, the only hydrogen line observed in the avail- 
able spectra. This line shows a well defined blue damping 
wing, but the red damping wing is partially blended with 
Hi clouds (see Fig. The Hi column density measure- 
ment was hampered by an additional difficulty, the location 
of the sub-DLA Lya line in the red wing of the DLA Lya 
line. Therefore, we first had to locally renormalize the quasar 
spectrum with the fit of the DLA Lya damping wing pro- 
file, before measuring the N(Hl) of the sub-DLA by fixing 
the redshift at the position of the strongest component of 
the metal line profiles (the component 6). The error on the 
derived hydrogen column density was estimated by varying 
the continuum level by 5%. 

As shown in Fig. [7| the low-ion transitions of the sub- 
DLA system present a complex velocity structure com- 
posed of 11 components spread over ~ 210 km s" 1 (see 
Table |2J . The numerous Fe + lines identified in the system 
(Fen A1608,2382, 2586,2600) allowed to obtain a very accu- 
rate A r (Fell) measurement. We also provide measurements 
of the Si 11 and Al 11 column densities. To be conservative, we 
prefer to consider the measured N(0 1) obtained from the 
Oi A1039,1302 lines located in the Lya forest as an upper 
limit. The Znn and Cm lines are not detected in the spec- 
tra, we provide 4 a upper limits. The C 11 lines are highly 
saturated. 

In addition, we observe the C IV and Si IV high- ion tran- 
sitions showing in this case very similar line profiles to the 
low-ion transitions. The C IV A1550 and Siiv A1402 lines are 
blended with the metal lines of other absorbers (see Fig. 01. 
We also detect the Aim and Fein intermediate- ion transi- 




Veloeity [km s ! ] 



Figure 8. Normalized UVES spectrum of Q1409+095 showing 
the sub-DLA Lya profile with the Voigt profile fit. The zero ve- 
locity is fixed at 2 = 2.668227. The vertical bar corresponds to the 
velocity centroid of the strongest component of the low-ion tran- 
sition lines used for the best fit, the component 2 at z = 2.668227. 
The measured Hi column density is logiV(Hl) = 19.75 ± 0.10. 

tions. The intermediate-ion column density determination, 
in particular the one of Fe ++ , is crucial for the estimation 
of the ionization corrections, which may not be negligible 
in this low Hi column density system (see Section [IJ. Since 
the Aim lines are weak and the Fein line is in the Lya 
forest where blends with Hi clouds are possible, we cannot 
determine with high accuracy whether the intermediate-ion 
profiles present differences with respect to the low-ion pro- 
files in this sub-DLA system. Thus, as the Aim lines show 
a single velocity component located exactly at the redshift 
of the strongest component of the low-ion line profiles, the 
component 6, we assume that we can adopt the same fitting 
parameters for the intermediate-ion transitions as the ones 
defined for the low-ion transitions. 

3.4 Q1409+095, z sub - D LA = 2.668 

Discovered by Hazard et al. (1986), this quasar has first been 
studied at low resolution by Lanzetta et al. (1991) and then 
by e.g. Lu et al. (1993). They reported the identification of 
a DLA system at z abs = 2.456 and another absorption line 
system at z abs = 2.668 with a Hi column density below the 
DLA definition corresponding to a sub-DLA system. Pettini 
et al. (2002) studied these two systems at high resolution 
and obtained the abundance measurements of O, Si, N and 
Fe. We identified the sub-DLA system in the UVES spectra 
thanks to the Lya absorption line with clearly visible damp- 
ing wings (see Fig. |SJ, and derived a Hi column density of 
log AT(Hl) = 19.75±0.10 in agreement with the Pettini et al. 
measurement. 

Our determinations of the Oi, Sill and Fell column 
densities and of the Ni upper limit are also consistent 
with their findings, although we obtained ~ 0.1 dex higher 
Sill and Fell column densities. This discrepancy may be 
explained by the fact that we used a two-component fit- 
ting model for the low-ion transitions (see Table 2J, while 
these authors used a one-component fitting model. Indeed, 
weak metal lines, like the Ni, Sn lines and some Fell lines 
with low oscillator strengths, are well fitted with a single 
component. However, as shown in Fig. [§] the strong Fell 
A2382,2586,2600, Oi A1302, Sin A1190 and Cn A1036,1334 
line profiles clearly indicate the presence of a second weaker 
component at ~ 11 km s _1 bluewards of the main compo- 
nent. In addition to the O , Si + , Fe + column densities and 
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Figure 9. Same as Fig. for the sub-DLA towards Q1409+095. The zero velocity is fixed at z = 2.668227. The vertical bars mark the 
positions of the fitted velocity components (see Table |41 . 



the N(N i) upper limit — the N° triplet being blended with 
the Lya forest absorptions — , we obtained the column den- 
sity measurements of S + from the Sn A1250,1259 lines and 
Cr+. The Znn lines are not detected, thus we provide a 4 a 
upper limit. We also report a lower limit on iV(Cll), since 
the feature is saturated. 

We detect the Al ++ and Fe ++ intermediate-ion tran- 
sitions in this sub-DLA system. A single component is ob- 
served in the Aim A1854 and A1862 line profiles, and it is 
reasonably well fitted with the low-ion transition fitting pa- 
rameters (see Fig. Although the fit of the Aim A1854 
component does not look optimal and could be shifted by a 
few km s _1 bluewards, when fitting simultaneously the two 
Aim lines at A rcst = 1854 and 1862 A, the best fit solution 
is very similar to the low-ion transition fitting parameters. 
Thus, given the weak intensity of the Al ill and Fe m absorp- 
tion lines and the low S/N in the spectral regions covering 
these lines, we preferred to adopt the same fitting param- 



eters as the ones derived for the low-ion transition lines, 
and we obtained Aim and marginal Fein column density 
measurements. No high-ion transition line is observed, all of 
them being outside the available wavelength coverage. 



3.5 Q1444+014, z sub -DLA = 2.087 

The quasar Q1444+014 is found more commonly under the 
name LBQS 1444+0126 and was discovered by Hewett et al. 
(1991). Its intervening absorption system at z a b s = 2.087 
was first reported by Wolfe et al. (1995). Pettini et al. (2002) 
analysed the H I, O I, Si II, N I and Fe II lines of this absorber, 
loosely referring to it as a DLA system. We derived the 
H I column density by fitting the weakly contaminated Lya 
absorption line showing well defined damping wings (the 
redshift was left as a free parameter; see Fig. 1101 . and we 
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Table 5. Component structure of the z a t, B = 2.087 sub-DLA system towards Q1444+014 



No z abs v* el 


K<r b ) 


Ion 


logAr(o- logJV ) 


No 


■^abs 


rcl 




Ion 


log iV(cr log N ) 


km s — 1 


km s _1 




cm 






km s — 1 


km s _1 




cm 


Low-ion transitions 


1 2.085060 -170 


10.9(1.0) 


Sin 


12.98(0.06) 


8 


2.086808 





13.3(0.2) 


Si II 


14.37(0.01) 






Oi 


13.75(0.04) 










Oi 


> 15.81 






Cn 


13.44(0.03) 










Cii 


> 15.64 


2 2.085297 -147 


5.2(0.6) 


Sin 


13.27(0.05) 










Aln 


12.66(0.01) 






Oi 


> 14.69 










Fen 


13.91(0.02) 






Cn 


> 13.70 


9 


2.087482 


+65 


7.1(0.2) 


Sin 


13.54(0.02) 






Alii 


11.88(0.08) 










Oi 


14.38(0.02) 






Fen 


12.81(0.08) 










Cii 


14.57(0.03) 


3 2.085481 -129 


10.4(0.5) 


Sin 


14.02(0.02) 










Aln 


12.27(0.02) 






Oi 


> 14.99 










Fen 


13.26(0.04) 






Cii 


> 16.07 


10 


2.087682 


+85 


6.3(0.8) 


Sin 


12.73(0.08) 






Alii 


12.81(0.02) 










Oi 


13.75(0.03) 






Fen 


13.66(0.02) 










C ii 


13.65(0.02) 


4 2.085711 -107 


6.1(0.4) 


Sin 


13.53(0.02) 










Aln 


11.07(0.19) 






Oi 


14.59(0.25) 


11 


2.087973 


+113 


5.6(0.4) 


Sin 


12.92(0.05) 






C ii 


> 14.26 










O I 


1 3 86co 041 






Alii 


12.06(0.05) 










Cii 


13.82(0.03) 






Fen 


13.28(0.03) 










Aln 


11.61(0.07) 


5 2.085857 -92 


8.5(3.2) 


Sin 


13.16(0.12) 


12 


2.088111 


+127 


6.1(0.4) 


Sin 


13.30(0.03) 






Oi 


14.20(0.12) 










Oi 


14.24(0.02) 






Cii 


> 13.95 










Cii 


14.32(0.21) 






Aln 


11.75(0.12) 










Aln 


12.01(0.04) 


6 2.086125 -66 


10.2(0.4) 


Sin 


14.09(0.03) 










Fen 


13.13(0.05) 






Oi 


> 15.36 


13 


2.088333 


+148 


13.8(2.0) 


Sin 


12.97(0.06) 






Cii 


> 15.24 










Oi 


13.75(0.05) 






Aln 


12.75(0.02) 










Cii 


14.03(0.10) 






Fen 


13.67(0.02) 










Aln 


11.75(0.07) 


7 2.086501 -30 


6.0(1.3) 


Sin 


12.81(0.07) 


11 


2.088617 


+176 


6.0(0.2) 


Sin 


13.66(0.02) 






Oi 


13.85(0.22) 










Oi 


14.69(0.04) 






Cii 


13.33(0.18) 










Cii 


> 14.95 






All! 


11.36(0.11) 










Aln 


12.48(0.02) 


















Fen 


13.40(0.02) 


High-ion transitions 


1 2.088197 +135 


17.5(1.8) 


Si iv 


12.69(0.04) 


2 


2.088584 


+173 


10.4(1.1) 


Si iv 


12.58(0.04) 






Civ 


13.05(0.04) 










Civ 


12.72(0.07) 


* Velocity relative to z = 2.086808 
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Figure 10. Normalized UVES spectrum of Q1444+014 showing 
the sub-DLA Lyo profile with the Voigt profile fit. The zero ve- 
locity is fixed at z = 2.086808. The vertical bar corresponds to 
the velocity centroid used for the best fit, z = 2.086643. The 
measured Hi column density is logiV(Hl) = 20.18 ± 0.10. 



classified the system as a sub-DLA with log ./V(Hl) = 20.18+ 
0.10, in agreement with the Pettini et al. measurement. 

The low-ion transition profiles of this sub-DLA system 



exhibit a very complex velocity structure with 14 compo- 
nents extended over ~ 350 km s" 1 (see Tableland Fig. ITU. 
We obtained column density measurements of Sill, Aln and 
Fell. Only the strongest components are observed in the Fell 
A1144,1608 line profiles. Unfortunately, several components 
between v ~ and —150 km s" 1 of the Oi A1302 and Cn 
A1334 lines are saturated, thus we could only provide lower 
limits on the column densities of these two metal transi- 
tions. No intermediate-ion line was detected. We therefore 
give only the 4 a upper limits on TV(Allll) and JV(Fem) — 
the Fein A1122 line is located in the Lya forest, but the 
spectral region where the strongest components of the Fe ill 
line should be detected is free from important Hi Mend- 
ings, and thus we can derive a reliable 4 a upper limit on 
./V(Felll). A very surprising particularity observed in this 
system is the behaviour of the high-ion transition lines, C iv 
and Si iv. They are characterized by only two relatively weak 
components and are extended over a much narrower velocity 
range of ~ 40 km s _1 than the low- ion profiles (see Table [K] 
and Fig. II IB . Most of the DLA systems show the opposite 
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Figure 11. Same as Fig.[5]for the sub-DLA towards Q1444+014. The zero velocity is fixed at z = 2.086808. The vertical bars mark the 
positions of the fitted velocity components (see Table 151 . 



situation with the high-ion line profiles spread over a larger 
velocity range than the low-ion line profiles. 



3.6 Q1451 + 123 

The absorption line systems towards this quasar discovered 
by Hazard et al. (1986) have first been studied by Wolfe 
et al. (1986). Lanzetta et al. (1991) identified in their low- 
resolution spectra a DLA system at z a b s = 2.469 and two 
other absorption systems at z a b s = 2.255 and 3.171. Pe- 
titjean, Srianand & Ledoux (2000) studied the intervening 
DLA and sub-DLA at z^hs = 3.171 for their molecular con- 
tent and report some of their elemental abundances. We 
present here for the first time the abundance analysis of the 
additional system at z a b s = 2.255 which is a borderline case 
between the DLA and sub-DLA systems. 



2 
E 




-3000 3000 

Velocity [km s -1 ] 

Figure 12. Normalized UVES spectrum of Q1451+123 showing 
the Lya profile with the Voigt profile fit of the borderline case 
between the DLA and sub-DLA systems. The zero velocity is fixed 
at z = 2.254710. The vertical bar corresponds to the velocity 
centroid used for the best fit, z = 2.255151. The measured Hi 
column density is logiV(Hl) = 20.30 ± 0.15. 



3.6.1 ZDLA/sub-DLA — 2.255 

This absorption metal-line system was identified thanks to 
the Lya absorption line observed only ~ 100 A away from 
the blue edge of the quasar spectrum (the quasar flux cut- 
off is at ~ 3835 A). The normalization is very difficult in 
this region, but the limited number of Lya forest absorp- 
tions leads to a clear high H I column density system detec- 
tion, the damping wings being well observed (see Fig. 1121 . 
We obtained a hydrogen column density measurement of 
log A r (H i) = 20.30 ± 0.15, which shows that this absorber is 



a borderline case between the DLA and sub-DLA systems. 
A large error bar on N{H i) was adopted to take into account 
the uncertainty in the continuum level. This absorption line 
system, although analysed in this paper, is not considered 
as a sub-DLA and is not included in the sub-DLA sample 
used for the different studies undertaken in Paper II. 

Table |S| and Fig. 1131 present the velocity structure of 
the low-ion transitions in this absorber. It is characterized 
by 5 components spread over ~ 130 km s _1 . We obtained the 
column density measurements of Zn + , Cr + and Fe + from the 
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Figure 13. Same as Fig. [F]for the borderline case between the DLA and sub-DLA systems towards Q1451+123. The zero velocity 
fixed at 2 = 2.254710. The vertical bars mark the positions of the fitted velocity components (see Table lt)l . 



Table 6. Component structure of the z a b s = 2.255 borderline 
case between the DLA and sub-DLA systems towards Q1451+123 



No z abs v* cl 


h(cr b ) 


Ion 


logAf(o- logJV ) 


km s _1 


km s _1 




cm -2 


Low-ion transitions 


1 2.254420 -27 


11.5(0.6) 


Sin 


> 14.59 






Fen 


13.68(0.02) 






Nin 


< 13.20 






Cm 


12.33(0.15) 


2 2.254584 -12 


5.0(2.1) 


Sin 


> 14.56 






Fen 


13.73(0.08) 






Nin 


< 12.77 






Znn 


11.49(0.15) 






Crn 


12.37(0.13) 


3 2.254710 


2.8(1.4) 


Sin 


> 15.31 






Fen 


13.84(0.12) 






Nin 


< 13.09 






Znn 


11.60(0.08) 






Crn 


12.48(0.15) 


4 2.255047 +31 


15.4(0.4) 


Sin 


13.95(0.08) 






Fen 


13.56(0.01) 


5 2.255725 +93 


7.2(0.7) 


Sin 


12.56(0.21) 






Fen 


12.86(0.03) 


Intermediate-ion transitions 


1 2.254395 -29 


18.2(3.5) 


Aim 


12.29(0.07) 


2 2.254675 -3 


5.1(0.6) 


Aim 


12.48(0.03) 



* Velocity relative to z = 2.254710 

Fen A2374,2382,2600 lines. Only the strongest components, 
the components 1, 2 and 3, are observed in the weak Znn 
and Cr II lines detected thanks to the relatively high H I 
column density of the system. Several interesting metal lines 
are unfortunately located in the spectral wavelength gap 
between 4855 and 5900 A. The only observed Si + line in the 



Table 7. Component structure of the z^ s = 3.171 sub-DLA 
system towards Q1451+123 



No 



v* el _ b(a b ) 
km s 1 km s 1 



Ion 



logiV((T logjv ) 
™-2 



Low-ion transitions 



1 3.170795 -23 


8.9(1.3) 


Sin 


13.29(0.04) 






Oi 


< 14.46 






Cn 


13.96(0.18) 






Aln 


11.97(0.03) 






Fen 


12.99(0.07) 


2 3.171109 


9.2(0.9) 


Sin 


13.38(0.03) 






Oi 


< 14.53 






Cn 


14.04(0.21) 






Aln 


12.04(0.03) 






Fen 


13.06(0.06) 


High-ion transitions 


1 3.172149 +75 


14.3(3.0) 


Civ 


13.39(0.14) 


2 3.172427 +95 


12.4(2.4) 


Civ 


12.78(0.07) 



Velocity relative to z = 3.171109 



available spectra is the saturated Sin A1304 line which is, in 
addition, blended with the damping red wing of the Lya line 
of the DLA system at z a b s = 2.469. To obtain a correct lower 
limit we first had to renormalize the portion of the quasar 
spectrum around the Sill line with the fit of the DLA Lya 
damping wing profile. The Nin A1317,1370 lines located in 
the Lya forest are detected at the noise level and at the 
uncertainty level of the continuum placement. We therefore 
considered the measured Ni + column density as an upper 
limit. The Oi A1302, Cn A1334 and the Mgn doublet lines 
are also detected, but these lines are so highly saturated that 
no useful lower limit can be deduced. 

From the higher ionization transitions, only the Al ++ 
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Table 8. Component structure of the z a b s 
system towards Q1511+090 



2.087 sub-DLA 



No 


■^abs 


v* , 

rcl 


b(a b ) 


Ion 


logAf(cr logiv ) 
cm -2 






km s _1 


km s _1 




Low-ion transitions 


1 


2.086977 


-121 


9.1(0.3) 


Sin 
Oi 
Alii 
Fen 


12.84(0.01) 
13.14(0.03) 

11.66(0.18) 

12.14(0.02) 


2 


2.087631 


-57 


30.0(3.8) 


Sin 
Oi 
Alii 
Fen 


12.87(0.02) 

< 13.81 
12.29(0.13) 
12.12(0.03) 


3 


2.087833 


-38 


11.4(1.0) 


Sin 
Oi 
Alii 
Fen 


12.70(0.03) 
< 13.06 

11.81(0.14) 

12.05(0.04) 


4 


2.088223 





13.5(0.8) 


Sin 
Oi 
Al ii 
Fen 


12.98(0.02) 

< 13.94 
11.78(0.14) 

12.19(0.03) 


5 


2.088643 


+41 


11.1(1.0) 


Sin 
Oi 
Alii 
Fen 


12.80(0.01) 

< 14.23 
11.76(0.16) 

12.12(0.03) 


6 


2.088861 


+62 


5.6(0.3) 


Sin 
Oi 
Alii 
Fen 


13.17(0.01) 

< 14.52 
11.99(0.13) 
12.37(0.02) 


7 


2.089054 


+81 


9.7(0.5) 


Sin 
Oi 
Alii 
Fen 


13.47(0.18) 

< 14.93 
12.24(0.13) 
12.73(0.03) 


8 


2.089310 


+106 


17.6(1.4) 


Sin 
Oi 
Alii 
Fen 


13.50(0.15) 

< 14.68 
12.22(0.13) 
12.58(0.04) 



* Velocity relative to z = 2.088223 

intermediate-ion lines and the Si +++ high-ion lines are cov- 
ered in the available spectra. The Aim lines clearly show 
slightly different profiles from the low-ion profiles, but with a 
very similar trend (see Fig. 1131 . We provide an Al ++ column 
density measurement obtained with the fitting parameters 
given in Table |S| No column density measurement can be 
obtained from the Siiv doublet, since it is heavily blended 
with Lya forest absorptions. 



3.6.2 Zsub-DLA = 3.171 

The Lya absorption line of this sub-DLA system is outside 
the available spectral coverage. We identified this sub-DLA 
system thanks to the Hi column density measurement re- 
ported by Petitjean, Srianand & Ledoux (2000). They esti- 
mated the hydrogen column density of this system from the 
equivalent width measurement by Bechtold (1994) and give 
log AT(H i) ~ 19.70. We confirm this value with the measure- 
ments obtained from the fits of the Ly/3, Ly7 and Ly8 lines 
(the redshift and the b- value were left as free parameters; see 
Fig. 1141 . The other members of the Lyman series were not 
taken into account, being heavily blended with Lya forest 
absorptions. 
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Figure 14. Normalized UVES spectrum of Q1451+123 showing 
the sub-DLA Ly/3, Ly7 and Ly8 profiles with the Voigt profile fits. 
The zero velocity is fixed at z = 3.171109. The vertical bar corre- 
sponds to the velocity centroid used for the best fit, z = 3.171606. 
The measured Hi column density is logAf(Hl) = 19.70 ± 0.15. 



This sub-DLA system shows a very simple velocity 
structure. The low-ion transition profiles are well fitted with 
two components (see Table|7|and Fig. 1151 . We obtained the 
column density measurements of Si + , C + , Al + and Fe + , and 
confirm the Petitjean, Srianand & Ledoux (2000) [Fe/H] and 
[Si/H] abundance measurements. The Cn column density 
was deduced from the C II A1036 line located in the Lya for- 
est. We adopted a relatively large error on TV(Cll) to take 
into account possible Hi contaminations. We observe eight 
O I lines in the available spectra, but the bulk of them are 
heavily blended with Lya forest absorptions. Only the O I 
A988,1039 lines allowed us to obtain a significant upper limit 
on iV(Ol). The Znn and Cm lines are not detected in the 
spectra, thus we provide only 4 a upper limits. No column 
density measurement of intermediate-ion transitions was ob- 
tained. The Fein A1122 line is heavily blended in the Lya 
forest, and the Aim lines are not detected, hence only a 4 a 
upper limit was derived. 

Similarly to the sub-DLA system at z a b s = 2.087 to- 
wards Q1444+014, the Civ high-ion transition lines ob- 
served in this sub-DLA show a peculiar behaviour. They 
are shifted by ~ 80 km s _1 redwards of the low-ion transi- 
tion lines. Although such a velocity difference between low- 
and high-ion lines is also observed in some DLA absorbers 
(Ledoux et al. 1998), in this particular case, the low- and 
the high-ion profiles do not overlap at all in velocity space 
(see Table [7| and Fig. [^J . This could suggest that the C iv 
lines are not associated with the studied sub-DLA system. 
Under this assumption, it may mean that another metal- 
line system at ~ +80 km s _1 is blended with the absorber 
associated with the low-ion transitions. However, this situ- 
ation seems unlikely, since we should then observe the ab- 
sorption lines of higher members of the Lyman series break- 
ing into two components. A weak emission feature is ob- 
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Figure 15. Same as Fig.[S]for the sub-DLA towards Q1451+123. The zero velocity is fixed at z = 3.171109. The vertical bars mark the 
positions of the fitted velocity components (see Table 171 . 



served at ~ 40 km s _1 in the Ly8 line (see Figure and it 
may cast some doubt on the previous statement. However, 
given the low quality of the data, it is difficult to determine 
with confidence whether this feature is real or just noise. 
Indeed, the feature is detected at less than 2 a in the Ly8 
line which is observed less than 15 A away from the blue 
edge of the quasar spectrum (the quasar flux cut-off is at 
~ 3835 A), where the signal-to-noise ratio is very low (only 
~ 5 FWHM). Therefore, with the available data we assume 
that the analysis of this absorption line system presented 
here is correct. However, we advise the reader that higher 
S/N spectra in the blue are required to determine whether 
the observed emission feature in the Ly8 line is real or not, 
since it will then imply a revision of the H I column density 
measurement in this absorber. If in the end this system is 
not a sub-DLA, it has no impact on our unbiased sample of 
sub-DLAs and on the statistical work presented in Paper II, 
since due to the non-observation of the Lya line in this sys- 
tem, it was excluded from the sample of sub-DLAs used for 
the statistical work (see Section \l. 31 . 



3.7 Q1511+090, z sub _ DLA = 2.088 

The absorption line systems towards this quasar were first 
analysed by Sargent, Steidel & Boksenberg (1988). But, the 
sub-DLA system at z a b s = 2.088 has never been previously 
studied, nor reported in the literature. It has been identified 
thanks to the Lya absorption line, the only hydrogen line 
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Figure 16. Normalized UVES spectrum of Q1511+090 showing 
the sub-DLA Lya profile with the Voigt profile fit. The zero ve- 
locity is fixed at z = 2.088223. The vertical bar corresponds to 
the velocity centroid used for the best fit, z = 2.088330. The 
measured Hi column density is log ./V(H i) = 19.47 ± 0.10. 



observed in the available spectra. The blue damping wing of 
the Lya line is blended with other Hi clouds, and only the 
outer parts of the blue wing are visible. The red damping 
wing is on the contrary relatively well constrained and can 
be used to determine the hydrogen column density of the 
system (see Fig. 11611 . By fixing the b- value at 20 km s _1 
and by leaving the redshift as a free parameter, we obtained 
a Hi column density of logAT(Hi) = 19.47 ± 0.10 as the 
best fit solution for the Lya line. The derived redshift at 
z = 2.08833 is located between the components 4 and 5 of 
the low-ion transition profiles. 
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Figure 17. Same as Fig. Elfor the sub-DLA towards Q1511+090. The zero velocity is fixed at z = 2.088223. The vertical bars mark the 
positions of the fitted velocity components (see Table 151 . 



Table |S| and Fig. 1171 present the velocity structure of 
this sub-DLA system. The low-ion transition lines are well 
described by 8 components extended over a large velocity 
interval of ~ 240 km s _1 . All the metal lines with rest- 
wavelengths up to Arcst = 1527 A are located in the Lya 
forest, thus the column density determinations of a large 
number of elements have to be done with caution. We ob- 
tained the column density measurements of Sill, Alii and 
Fen. The access to the strong Fen A2344,2382,2600 transi- 
tions allowed us to measure accurately A r (Fe n) of this rather 
low metallicity system, [Fe/H] = -1.71±0.10. The O I A1302 
line is blended with Lya forest lines, thus no reliable column 
density measurement can be determined, and we choose to 
consider the derived value as an upper limit. The case of the 
Cn A1334 line is even more difficult. The saturation of the 
line together with the strong blending, makes impossible the 
estimation of the column density of C + . The Al ++ and Fe ++ 
intermediate-ion transition lines are not detected, hence we 
provide only 4 a upper limits on A^(Allll) and A r (Fem) — 
the Fein A1122 line is located in the Lya forest, but the 
spectral region where the strongest components of the Fe ill 
line should be detected is free from Hi blendings, and thus 
we can derive a reliable 4 a upper limit on A^(Fe ill). The ob- 
served high-ion transition lines, the C IV and Si IV doublets, 
are all heavily blended with Lya forest lines preventing us 
from measuring their column densities. 

3.8 Q2059-360, z sub -DLA = 2.507 

The quasar Q1451+123 was discovered by Warren, Hewett 
& Osmer (1991), and it has been studied for its intervening 
DLA system at z a b s = 3.083. The molecular content of this 



Table 9. Component structure of the z a |-, a 
system towards Q2059— 360 



2.507 sub-DLA 



No 2 abs v* el 


b(<r b ) 


Ion 


log Af(o-i og jv) 


km s _1 


km s — 1 




cm 


Low-ion transitions 


1 2.507341 


4.9(0.2) 


Sill 


13.66(0.08) 






Oi 


15.52(0.21) 






C ii 


14.82(0.21) 






N i 


12.66(0.18) 






Fen 


13.41(0.03) 


2 2.507535 +17 


5.3(0.8) 


Sin 


12.91(0.25) 






Oi 


14.11(0.13) 






Cn 


13.87(0.15) 






Feu 


12.59(0.04) 



* Velocity relative to z = 2.507341 

DLA has been analysed by Petitjean, Srianand & Ledoux 
(2000). We have found in the spectra of Q2059-360 an addi- 
tional high hydrogen column density absorption line system, 
a sub-DLA at z a b s = 2.507. It has been identified thanks to 
the Lya absorption line which shows well defined damping 
wings (see Fig. 1181 . No other hydrogen line is observed in 
the available spectra. The fit of the Lya line was performed 
by fixing the fe- value at 20 km s _1 and the redshift at the 
redshift of the strongest component of the low-ion transition 
lines, the component 1. The derived Hi column density is 
log AT (Hi) = 20.21 ±0.10. 

The low-ion transition profiles are characterized by a 
simple velocity structure composed of two components (see 
Table |5| and Fie. 1191 . In the weaker lines, only the strongest 
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Figure 19. Same as Fig.0for the sub-DLA towards Q2059-360. The zero velocity is fixed at z = 2.507341. The vertical bars mark the 
positions of the fitted velocity components (see Table 151 . 
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Figure 18. Normalized UVES spectrum of Q2059-360 showing 
the sub-DLA Lya profile with the Voigt profile fit. The zero veloc- 
ity is fixed at z = 2.507341. The vertical bar corresponds to the 
velocity centroid of the strongest component of the low-ion tran- 
sition lines used for the best fit, the component 1 at z = 2.507341. 
The measured Hi column density is logAT(Hl) = 20.21 ± 0.10. 



component at v — is observed. We obtained an accurate 
Fe + column density measurement thanks to the numerous 
observed Fell lines at A rcst = 2344, 2374, 2382, 2586 and 
2600, and despite the fact that the reddest parts of the avail- 
able spectra at A > 8000 A are noisy with a S/N per pixel 
lower than 10. The observed Sill A1304, Ol A1302 and Cn 



A 1334 lines are all located in the Lya forest. Large errors 
have been adopted on their column density measurements to 
take into account possible Lya contaminations. In the case of 
Si + , we could constrain the column density of the strongest 
component with the Sill A1808 line located outside the Lya 
forest. However, the derived O and C + column density 
measurements might be considered as borderline cases be- 
tween values and limits. We also obtained a column density 
measurement of N° from the Nl A1200. 223,1200. 710 lines, 
the Nl A1199.550 line being blended. This measurement is 
marginal, since the detected N I lines are extremely weak, at 
the limit of the noise level. The Zn II and Cr II lines are not 
detected despite the relatively high Hi column density of 
the system, thus we provide 4 a upper limits. We also give 
the 4 a upper limit on the Al ++ intermediate-ion transition 
column density, the Aim lines being not detected. No high- 
ion transition column density has been measured, since the 
Si IV lines are heavily blended with Lya forest absorptions 
and the C IV lines are not covered in the available spectra. 

3.9 Q2116-358, z 8ub -DLA = 1-996 

The quasar Q2116— 358 was discovered by Adam (1985). Its 
intervening absorption metal-line system at z a b s = 1.996 
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Figure 21. Same as Fig.0for the sub-DLA towards Q2116-358. The zero velocity is fixed at z = 1.996119. The vertical bars mark the 
positions of the fitted velocity components (see Table [TU1 . 
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Figure 20. Normalized UVES spectrum of Q2116— 358 showing 
the sub-DLA Lyo profile with the Voigt profile fit. The zero ve- 
locity is fixed at z = 1.996119. The vertical bars correspond from 
right to left to the velocity centroids of the component 6 of the 
low-ion transition lines at z = 1.996224 and z = 1.992317 used 
for the best fit, belonging to the sub-DLA system and an addi- 
tional absorber, respectively. The dotted line corresponds to the 
fit with the sub-DLA only. The measured H I column density of 
the sub-DLA is logTV(Hl) = 20.06 ± 0.10. 



has previously been reported by Lanzetta, Wolfe & Turn- 
shek (1987), Wampler, Bergeron & Petitjean (1993), M0ller, 
Jakobsen & Perryman (1994), and it has even been targeted 
in an imaging programme by Warren et al. (2001). In this 
paper we provide a new complete abundance analysis of this 
system using the UVES high-resolution spectra. 

Only the hydrogen Lya line is observed in the available 
data. It shows well defined damping wings which allowed an 
accurate measurement of the H I column density of the sys- 
tem. We obtained a value of log N (Hi) = 20.06 ±0.10, lower 
by ~ 0.1 dex than the value reported by M0ller, Jakobsen 
& Perryman (1994), and which indicates that the absorber 
is a sub-DLA system. Fig. 1201 shows the best fit solution of 
the Lya line. The presence of an asymmetry in the profiles 
of the blue and red damping wings made it necessary to 
include, besides the contribution of the sub-DLA system, a 
second low column density absorption system bluewards of 
the sub-DLA in order to obtain a good fit of both wings 
of the Lya line. During the fitting procedure, the redshift 
and the 6-value of the sub-DLA system were fixed at the 
redshift of the component 6 of the low-ion transition lines 
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Table 10. Component structure of the 2 a b s = 1-996 sub-DLA system towards Q2116— 358 



No z abs v* ei 


&fo) 
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log JV(CTlogjv) 
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Ion 


log N(cr logN ) 


km s _1 


km s _1 




cm~ 2 






km s 


km s _1 




cm -2 


Low-ion transitions 


1 1.995416 -70 


5.0(0.8) 


Si II 


14.39(0.07) 


6 


1.996224 


+10 


7.5(2.6) 


Si ii 


14.71(0.14) 






Oi 


> 15.69 










Oi 


> 16.00 






Sn 


< 14.12 










Sii 


< 14.42 
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Nin 


12.79(0.07) 
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Oi 
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Oi 
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Fern 


< 13.49 


2 1.995903 -22 


17.1(2.9) 


Aim 


12.67(o.io) 


5 


1.996458 


+34 


6.4(2.5) 


Aim 


11.87(0.11) 






Fein 


< 14.06 










Fein 


< 13.62 


3 1.996172 +5 


13.5(3.0) 


Aim 


12.80(o.ii) 


6 


1.996814 


+69 


5.5(0.6) 


Aim 


12.40(0.02) 






Fein 


< 13.99 










Fein 


< 13.51 


High-ion transitions 


1 1.994594 -153 


12.5(0.7) 


Si iv 


12.88(0.02) 


6... 


1.996344 


+22 


9.7(0.8) 


Si iv 


13.53(0.03) 






Civ 


13.33(0.06) 










Civ 


13.89(0.14) 


2 1.994974 -115 


16.3(i.i) 


Si iv 


12.90(0.02) 


7... 


1.996517 


+40 


6.4(0.8) 


Si iv 


12.65(0.05) 






Civ 


13.41(0.04) 










Civ 


13.48(0.16) 


3 1.995426 -69 


13.6(0.3) 


Si iv 


13.29(o.oi) 


8... 


1.996789 


+67 


13.7(0.5) 


Si iv 


13.81(0.01) 






Civ 


13.80(0.02) 










Civ 


> 15.03 


4 1.995957 -16 


17.4(1.3) 


Si iv 


13.41(0.05) 


9... 


1.997112 


+99 


13.2(0.5) 


Si iv 


13.44(0.01) 






Civ 


13.55(0.07) 










Civ 


> 14.43 


5 1.996150 +3 


9.1(0.6) 


Si iv 


13.52(0.04) 


10... 


1.997414 


+130 


9.7(0.4) 


Si iv 


13.22(0.02) 






Civ 


13.92(0.13) 










Civ 


14.33(0.12) 



* Velocity relative to z = 1.996119 
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and at 20 km s _1 , respectively, and the redshift and the 
b- value of the second absorber were left as free parameters. 
We obtained the best fit solution with the second absorption 
line system at z = 1.992317 shifted by about —391 km s _1 
from the sub-DLA system and with a Hi column density of 
logJV(Hi) = 14.59 ± 0.10. 

Metal lines from a wide range of ionization species are 
present in this sub-DLA system. Fig. 1211 shows the low-, 
intermediate- and high-ion transition line profiles, and their 
velocity structures are described in Table EU The low -ion 
transition lines are characterized by a relatively complex ve- 
locity structure of 11 components spread over ~ 220 km s . 
Only the strongest components, the components 3 to 6, are 
observed in the weak metal lines, such as the Nill, Cm and 
Znll lines. We obtained the column density measurements 
of Si + , Ni + , Cr + , Zn + and Fe + . Many components of the O I 
A1302 and Alii A1670 lines are strongly saturated, while the 
weaker components are well matched by the model param- 
eters defined from the Sill and Fell lines. Since no reliable 
value of the total column density of O and Al + can be de- 
duced, we provide only lower limits. We detect also the S II 
A1250,1259 lines. But, since they are blended with Lya for- 
est absorptions, we preferred to consider the measured S + 
column density as an upper limit. 

Relatively strong Aim lines are observed in this sub- 
DLA system (see Fig. 1211 . These intermediate-ion transition 
lines show profiles which are very similar to the low-ion pro- 
files, but need to be fitted with slightly different model pa- 
rameters given in Table tTol The Fe ill A1122 intermediate-ion 
line is also observed. Nevertheless, it is located in the Lya 
forest and is blended, therefore we provide only an upper 
limit on its column density. The observed high-ion transition 
lines are very strong and show a complex velocity structure 
composed of 10 components. They extend over a wider ve- 
locity range of ~ 290 km s _1 than the low-ion profiles (see 
Fig. I2H . We obtained the column density measurement of 
Si IV and a lower limit on N(G iv) — two components of the 
C IV doublet are heavily saturated. 



3.10 PSS J2155+1358 

High signal-to-noise ratio, 5 A resolution (FWHM) spectra 
of this Zem > 4 quasar have been obtained by Peroux et al. 
(2001). They identified a high Hi column density DLA sys- 
tem at z a b s = 3.316. The metallicity of this DLA system was 
first reported by Dessauges-Zavadsky et al. (2001b). Here we 
present the chemical analysis of three newly identified sub- 
DLA systems at z a b s = 3.142, 3.565 and 4.212 towards this 
quasar. 



3.10.1 z sub -DLA = 3.142 

This sub-DLA system has been identified thanks to the Lya 
absorption line. No other hydrogen lines are observed in 
this sub-DLA, since they are beyond the quasar flux cut- 
off. Its Hi column density is thus poorly constrained, all 
the more since the Lya line lies in the far blue end of the 
quasar spectrum. Due to the high number of H I absorptions 
present in these spectral regions, the Lya damping wings 
are strongly blended and the continuum placement is uncer- 
tain. Fig. 1221 shows the best fit solution obtained by fixing 
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Figure 22. Normalized UVES spectrum of PSS J2155+1358 
showing the sub-DLA Lya profile with the Voigt profile fit. The 
zero velocity is fixed at z = 3.141998. The vertical bar corre- 
sponds to the velocity centroid of the strongest component of the 
low-ion transition lines used for the best fit, the component 1 
at z = 3.141998. The measured Hi column density is logTV(Hl) 
= 19.94 ±0.10. 



Table 


11. Component structure of the 




3.142 sub-DLA 


system towards PSS J2155+1358 






No 


2abs V* el b(u b ) 


Ion 


logAT(o- logJV ) 




km s _1 km s — 1 




cm -2 


Low- 


and intermediate-ion transitions 






1 


3.141998 5.7(1.0) 


Sin 


13.36(0.06) 






Oi 


< 14.22 






Cn 


14.05(0.20) 






Alii 


12.07(0.07) 






Fen 


12.77(0.20) 






Aim 


< 11.57 


2 


3.142313 +23 14.5(2.6) 


Sin 


12.76(0.10) 






Oi 


< 13.91 






Cn 


< 13.66 






Aln 


11.50(0.18) 


3 


3.142712 +52 10.0(2.8) 


Sin 


13.20(0.07) 






Oi 


< 14.42 






Cn 


< 14.00 






Aln 


11.81(0.13) 






Fen 


13.04(0.18) 



* Velocity relative to z = 3.141998 



the b- value at 20 km s^ 1 and the redshift at the redshift of 
the strongest component, the component 1, of the low-ion 
transition lines. We derived a hydrogen column density of 
log AT (Hi) = 19.94 + 0.10. 

Few metal lines were detected in this sub-DLA system, 
mainly because all of the lines with rest-wavelengths less 
than 1543 A are located in the Lya forest. The velocity 
structure of the low-ion transition lines was deduced from 
the Sill, Ol, Cn, Aln and Fell lines, all observed in the 
Lya forest apart from the Aln and Fell lines. It is well 
described by three components (see Table fTTI and Fig. 1231 . 
We were able to measure the column density of Si + thanks 
to the detection of three Sill lines at A rcs t = 1190, 1193 
and 1304 A. We obtained the column density measurement 
of Al + from the Aln A1670 line whose velocity components 
are just not blended with telluric lines, which we identified 
with the help of the spectrum of a hot, fast rotating star 
observed the same night as the quasar. We deduced the Fe + 
column density from the Fell A1608 line. The weakness of 
this Fe II line observed in this case just above the noise level 
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Figure 23. Same as Fig. for the sub-DLA towards PSS J2155+1358. The zero velocity is fixed at z = 3.141998. The vertical bars 
mark the positions of the fitted velocity components (see Table HT1 . 



and which is usually almost saturated in the DLA systems 
indicates that the metallicity in this sub-DLA is low, [Fe/H] 
= —2.21 ±0.21. To be conservative, we preferred to consider 
the derived O and C + column density measurements as 
upper limits due to the possible blends with Lya forest lines. 
No intermediate-ion and high-ion transitions are observed 
in this sub-DLA system. We provide an upper limit on the 
Al ++ column density obtained from the fit of the very weak 
feature observed in the Aim A1854 line at the redshift of 
the strongest component of the low-ion line profiles. 

3.10.2 Zs U b-DLA = 3.565 

This sub-DLA system was also identified thanks to the Lya 
absorption line. No other hydrogen line is observed in the 
available spectra. Given the high redshift of the quasar PSS 
J2155+1358, the Lya forest shows many absorptions, and 
thus it is not surprising that the damping wings of the Lya 
line are strongly blended. We made several fitting tests to 
find the best fit solution of the Lya line. The first test was 
done by fixing in the fitting procedure the 6-value at 20 
km s _1 and by leaving the redshift as a free parameter. We 
obtained a good fit with logAf(Hl) = 19.62 ± 0.15, but it 
was not a satisfactory solution, since the deduced redshift 
was shifted by about 73 km s _1 bluewards of the redshift 
of the single component observed in the low-ion transition 
lines. The second and adopted test was made by fixing both 
the b- value at 20 km s _1 and the redshift at the redshift 
of the single component of the low-ion lines. We obtained 
the fit shown in Fig. 1241 But, in this case, to fill the whole 
observed Lya absorption profile, we had to include, besides 
the sub-DLA contribution, a second low H I column density 
absorption system bluewards of the sub-DLA system. In the 
fitting procedure, the b- value and the redshift of this second 
absorber were left as free parameters. We obtained the best 
fit solution with the second absorption line system located 
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Figure 24. Normalized UVES spectrum of PSS J2155+1358 
showing the sub-DLA Lya profile with the Voigt profile fit. The 
zero velocity is fixed at z = 3.565272. The vertical bars correspond 
from right to left to the velocity centroids of the single component 
of the low-ion transition lines at z = 3.56527 and z = 3.561814 
used for the best fit, belonging to the sub-DLA system and an 
additional absorber, respectively. The dotted line corresponds to 
the fit with the sub-DLA only. The measured H I column density 
of the sub-DLA is logiV(Hi) = 19.37 ± 0.15. 



at z = 3.561814, i.e. shifted by about —228 km s _1 from the 
sub-DLA system, and with a H I column density of log A r (H i) 
= 16.46 ±0.15. The derived hydrogen column density of the 
sub-DLA is log A?"(H i) = 19.37 ± 0.15. 

As mentioned above, a very simple velocity structure 
with a single component is observed in the low-ion line pro- 
files (see Fig. 1251 and Table IT2^I . We obtained the column 
density measurements of Si + , C + and Al + . A r (Sill) was de- 
duced from four Sill lines at A rcst = 1190, 1260, 1304, 1526 
A. N(C n) is a marginal measurement, since the Cn A1334 
line is located in the Lya forest, where blends are not ex- 
cluded. Finally, Af(Aln) was obtained from the Alii A1670 
line found in a spectral region heavily contaminated by tel- 
luric lines. The spectrum of a hot, fast rotating star observed 
the same night as the quasar allowed us to identify the tel- 
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Table 12. Component structure of the 2 a bs — 

3.565 sub-DLA 

system towards PSS J2155+1358 



No 2 abs v* el 


b(<r b ) 


Ion 


logAf(o- logJV ) 


km s _ 


1 km s _1 




cm -2 


Low- and intermediate-ion transitions 


1 3.565272 


5.9(0.9) 


Sin 


13.66(0.04) 






Cn 


14.34(0.12) 






Alii 


12.16(0.12) 






Fc in 


< 13.81 






Aim 


12.71(0.09) 


High-ion transitions 


1 3.564519 -49 


12.5(2.6) 


Si iv 


12.95(0.07) 






Civ 


13.54(0.07) 


2 3.564939 -22 


10.8(i.i) 


Si iv 


13.38(0.03) 






Civ 


13.83(0.17) 


3 3.565289 +1 


7.4(1.2) 


Si iv 


13.57(0.04) 






Civ 


13.84(0.06) 


4 3.565531 +17 


5.5(0.9) 


Si iv 


13.02(0.04) 






Civ 


13.15(0.10) 


5 3.565859 +39 


3.7(i.i) 


Si iv 


12.97(0.08) 






Civ 


13.01(0.15) 


6 3.565966 +46 


18.8(1.0) 


Si iv 


13.36(0.03) 






Civ 


13.64(0.04) 


* Velocity relative to z = 


3.565272 






luric lines and to show that the Al + line is 


not blended, but 



the uncertainty relative to the continuum placement in this 
region remains. We therefore adopted a relatively large error 
on the measured AT(Aln). Only the Fell A1608 line is avail- 
able to measure the Fe + column density, and no absorption 
feature at the redshift of the single component observed in 
the other low-ion transition lines is detected. We thus pro- 
vide only a 4 (j upper limit on the Fe + column density, which 
indicates that the metallicity is very low in this sub-DLA 
system, lower than [Fe/H] < —2.31. 

The column density measurements of the intermediate- 
ion transitions are crucial to constrain the ionization param- 
eters in this sub-DLA system which has among the lowest 
Hi column density in our sample. We detect only the Aim 
A1862 line. The Aim A1854 line is blended with emission 
lines from badly subtracted cosmic rays in the reduction pro- 
cedure — only two exposures were obtained for this quasar. 
The observed Aim A1862 line shows a single component at 
exactly the redshift of the component of the low-ion tran- 
sition lines (see Fig. 1251 . Assuming that this feature is real 
and is Al ++ , we obtained a column density measurement 
using the same fitting parameters as those deduced for the 
low-ion lines. We also derived an upper limit on N(Fe ++ ) 
from the Fein A1122 line located in the Lya forest. In addi- 
tion, we detect the strong C iv and Si iv high- ion transition 
lines. They show a much more complex velocity structure, 
composed of 6 components and spread over ~ 120 km s _1 
than the low-ion lines (see Fig.^and Table [HI. 

3.10.3 ^ U b-DLA = 4.212 

This high redshift sub-DLA system was easily identified with 
the Lya absorption line which shows a very well defined red 
damping wing, the blue wing being partially blended with 



Table 13. Component structure of the z a b s = 4.212 sub-DLA 
system towards PSS J2155+1358 



No z abs v* el 


H<r b ) 


Ion 


logiV((T logJ v) 


km s - 


1 km s _1 




_9 

cm 


Low-ion transitions 


1 4.212229 


6.8(0.5) 


Sin 


13.01(0.04) 






Oi 


14.26(0.05) 






Cn 


13.70(0.07) 






Alii 


11.62(0.16) 






Fen 


12.76(0.20) 


2 4.212628 +23 


6.3(0.6) 


Sin 


12.88(0.04) 






Oi 


14.13(0.05) 






Cn 


13.59(0.05) 






Alii 


11.62(0.17) 






Fen 


12.43(0.30) 


High-ion transitions 


1 4.212394 +10 


22.5(6.7) 


Si iv 


12.87(0.10) 


2 4.212925 +40 


7.9(4.9) 


Si iv 


12.05(0.25) 


* Velocity relative to z = 


4.212229 








-1000 1000 

Velocity [km s -1 ] 



Figure 26. Normalized UVES spectrum of PSS J2155+1358 
showing the sub-DLA Lya, Lye and Lyll profiles with the Voigt 
profile fits. The zero velocity is fixed at z = 4.212229. The ver- 
tical bar corresponds to the velocity centroid used for the best 
fit, z = 4.212442. The measured Hi column density is logiV(Hl) 
= 19.61 + 0.10. 

Lya forest lines. The Lyman series is observed down to Lyl3 
in this sub-DLA and allowed us to obtain a very accurate 
Hi column density measurement by fitting simultaneously 
the Lya line and the unblended Lye and Lyll lines (see 
Fig. 1261 . In the fitting procedure, we left the redshift as a 
free parameter and we fixed the b- value at 23 km s _1 . The 
best fit solution was obtained with z = 4.212442 which fell 
just between the two components observed in the low-ion 
transition lines, and with a H I column density of log N(R i) 
= 19.61 ± 0.10. However, it can be seen that the sub-DLA 
system is not sufficient to fill the whole observed Lya ab- 
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Figure 25. Same as Fig. for the sub-DLA towards PSS J2155+1358. The zero velocity is fixed at z = 3.565272. The vertical bars 
mark the positions of the fitted velocity components (see Table If2l . 
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Figure 27. Same as Fig. for the sub-DLA towards PSS J2f 55+f 358. The zero velocity is fixed at z = 4.2f2229. The vertical bars 
mark the positions of the fitted velocity components (see Table [T3l . 
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Table 14. Component structure of the 2 a bs — 

3.882 sub-DLA 

system towards PSS J2344+0342 



No 


■^abs 


v* , 

rcl 


b(a b ) 


Ion 


logAT(o-i ogJV ) 






km s" 1 


km s _1 




cm" 2 


Low-ion transitions 


1 


3.880511 


-131 


8.7(2.0) 


Sin 


12.91(0.11) 










Cn 


13.55(0.06) 










Alii 


11.28(0.15) 


2 


3.880811 


-112 


5.6(1.3) 


Sin 


12.87(0.10) 










Cn 


13.48(0.06) 










Alii 


11.68(0.25) 


3 


3 881785 


—52 


7 4M 1 1 


Sin 


1 3 08m 051 










Cn 


13.79(0.06) 










Al ii 


11.92(0.09) 


4 


3.882032 


-37 


6 5(1 9) 


Si n 


1 9 fi^cn n7~\ 










P TT 


10.0U(U. 14) 


5 


3.882392 


-15 


4.2(3 3) 


C ii 




5 


3.882639 





6 4n 1 i 


Sin 


1 2 95ro 101 










C ii 












Al n 


1 1 ^Qtn i ^ 


7 


3.883733 


+67 




C ii 




g 


3.884153 


+93 


7 1 41 
t .1(3.4) 


C ii 




9 


3.884580 


+119 


10.5(0.9) 


Sin 


13.45(0.04) 










Cn 


14.12(0.04) 










Aln 


11.92(0.09) 


10 


3.885067 


+149 


12.1(1.2) 


Sin 


13.35(0.05) 










Cn 


14.09(0.03) 










Aln 


12.16(0.06) 


11 


3.885412 


+170 


3.5(2.0) 


Cn 


13.02(0.09) 


Hig 


h-ion transitions 








1 


3.883235 


+37 


20.0(5.7) 


Civ 


13.15(o.io) 


2 


3.883706 


+66 


11.8(2.8) 


Si iv 


12.91(0.07) 










Civ 


13.50(0.04) 


3 


3.884165 


+94 


10.7(1.4) 


Si iv 


13.00(0.04) 










Civ 


13.40(0.05) 


4 


3.884840 


+135 


18.8(3.9) 


Si iv 


13.28(o.io) 










Civ 


13.79(0.04) 


5 


3.885314 


+164 


7.7(1.4) 


Si iv 


12.81(0.06) 










Civ 


13.35(0.06) 



* Velocity relative to z = 3.882639 
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Figure 28. Normalized UVES spectrum of PSS J2344+0342 
showing the sub-DLA Lya profile with the Voigt profile fit. The 
zero velocity is fixed at z = 3.882639. The vertical bars corre- 
spond from right to left to the velocity centroids z = 3.884697 
and 2 = 3.881604 used for the best fit, belonging to the sub-DLA 
system and an additional absorber, respectively. The dotted line 
corresponds to the fit with the sub-DLA only. The measured H I 
column density of the sub-DLA is logJV(Hl) = 19.50 ± 0.10. 

the metallicity in this sub-DLA system is very low, [Fe/H] 
= -2.17 + 0.25. 

The Aim intermediate-ion transition lines are covered 
in the available spectra, but the signal-to-noise ratio is so 
low in this particular region that no Al ++ line is detected, 
and no reliable 4 a upper limit can be derived. In addi- 
tion, the Fein A1122, Si in A1206 and Cm A977 lines are 
all blended with Lya forest absorptions. Consequently, no 
intermediate-ion column density measurement was obtained 
in this sub-DLA system. On the other hand, the Si iv high- 
ion transition lines are observed (the C IV lines are outside 
the available spectral coverage). Again they show a different 
velocity structure from the low-ion line profiles, but are in 
this case extended over a similar velocity range (see Fig. \H\ 
and Table [Toll. 



sorption, so contributions from additional Lya clouds have 
to be included. 

Given the high signal-to-noise ratio of the available 
spectra in the regions where we detect the low-ion transition 
lines and the fact that they show a simple velocity structure 
composed of two well separated components (see Fig. 1271 
and Table Hoi, we obtained very accurate column density 
measurements of O , C + , Si + and Al + . Two Fell lines at 

= 1144 and 1608 A are observed in the available spec- 
tra. A very weak feature is detected in the Fe II A1144 line at 
the redshift of the component 1. But, the second component 
is partially blended with an emission line from badly sub- 
tracted cosmic rays in the reduction procedure — a single 
exposure was obtained for this quasar in the blue. No ab- 
sorption feature is observed in the Fe II A1608 lines, moreover 
the second component is blended with a badly subtracted 
strong sky line. We thus obtained only a marginal measure- 
ment of the Fe + column density, on which we adopted a 
large error. The weakness of these Fell lines indicates that 



3.11 PSS J2344+0342, z sub -DLA = 3.882 

Peroux et al. (2001) obtained 5 A resolution (FWHM) spec- 
tra of this quasar and identified a DLA system at z a b s = 
3.220. The metallicity of this DLA system was first reported 
by Dessauges-Zavadsky et al. (2001b). Here we present the 
chemical analysis of a newly discovered high-redshift sub- 
DLA system at 2 a b s = 3.882. It has been identified thanks 
to the Lya absorption line which shows partially blended 
damping wings. The Ly/3 and Ly7 lines are also observed 
in this sub-DLA system. The simultaneous fit of both Lya 
and Ly/3 — the Ly7 line is strongly blended — allowed us to 
obtain an accurate Hi column density measurement. 

A first series of fitting tests of the Lya and Ly/3 lines 
was done by considering a single absorption line system. By 
leaving the redshift and the fo-value as free parameters, the 
best fit solution was obtained with a sub-DLA system at 
z — 3.884697 — very close to the redshift of the strongest 
component of the low-ion transition lines, the component 
9 - and a Hi column density of log jV(Hi) = 19.50 ± 0.10. 
In this solution the red wings of the Lya and Ly/3 lines were 
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Figure 29. Same as Fig. 0for the sub-DLA towards PSS J2344+0342. The zero velocity is fixed at z = 3.882639. The vertical bars 
mark the positions of the fitted velocity components (see Table Hit . 



well fitted, but the fit clearly shows that another strong ab- 
sorption line system is present bluewards of the sub-DLA 
system (see Fig. 1281 . We then made a second series of fit- 
ting tests by considering this time two absorption line sys- 
tems. By leaving in the fitting procedure all the b- values and 
redshifts as free parameters, we obtained the best fit solu- 
tion with a sub-DLA system characterized by the values de- 
scribed above plus a second strong absorber at z = 3.881604 
shifted by about —190 km s _1 from the sub-DLA system 
with a Hi column density of log AT(Hl) = 18.94±0.15. This 
absorption line system is a borderline case between the Ly- 
man limit and sub-DLA systems, but we do not include it 
in the sub-DLA sample. 

Table 1141 and Fig. 1291 present the velocity structure of 
the sub-DLA system. The low-ion transition lines show a 
complex velocity structure well described by 11 components. 
We obtained accurate column density measurements of C + 
and Si + . The Alii A1670 line is located in a region strongly 
contaminated by telluric lines. We identified the telluric lines 
with the help of the spectrum of a hot, fast rotating star ob- 
served the same night as the quasar, and noticed that the 
components 2 and 7 are blended. We adopted a large error 
on the deduced AT(A1 + ) to take into account the contribu- 
tions of these blends. Two Fe + lines, Fell A1144 and Fell 
A1608, are observed. However, the Fell A1144 line is partially 
blended with Lya forest lines, and none of the 11 compo- 
nents is detected in the Fen A1608 line. We thus provide a 
4 a upper limit on the Fe + column density which indicates 
that the metallicity in this sub-DLA system is low, lower 
than [Fe/H] < —1.98. No intermediate-ion transition line is 
clearly observed. The Fein A1122, Cm A977, Nil A1083 and 
Sim A1206 lines are strongly blended with Lya forest lines, 
and the Al ill lines are not detected. We provide a 4 a upper 
limit on the Al ++ column density. 

It should be emphasized that the observed low-ion tran- 
sition line profiles extend over a very large velocity interval 



of ~ 320 km s , from -140 to +180 km s _1 in velocity 
space. It implies that the low-ion transition lines overlap 
over the two identified clustered absorption line systems — 
the sub-DLA and the absorber at z — 3.881604 shifted by 
only ~ —190 km s" 1 from the sub-DLA (see above). It may 
thus be possible that some of the components of the detected 
low-ion metal lines are associated with this second absorber 
and not with the sub-DLA system. Two distinct clumps in 
the low-ion profiles are not yet observed, but the two high 
H I column density absorption line systems are so close that 
it would not be surprising to find their metal lines blended. 
Moreover, the observed high-ion transition lines, C IV and 
Si IV, are surprisingly confined to a much narrower velocity 
interval from +25 to +170 km s" 1 (see Fig.|H]and Tablet!!!. 
This may suggest that the detected high-ion lines are asso- 
ciated only with the sub-DLA system. If this is the case, the 
components of the low-ion transition lines belonging to the 
sub-DLA could be only the ones which are located in the 
same velocity interval as the high-ion lines, i.e. the compo- 
nents 7 to 11. The components 1 to 6 could then be associ- 
ated with the second absorption line system. It is difficult to 
be categorical on the real situation, therefore it would per- 
haps be more sensible to speak of an interval of values for 
the column density measurements in this sub-DLA system. 



4 IONIZATION CORRECTION STUDY 

In the 12 sub-DLA systems and 1 borderline case between 
the DLA and sub-DLA systems studied, we obtained ionic 
column density measurements of low-ion transitions, O , 
C + , Si+, N°, S+, Mg+, A1+, Fe+, Ni+, Zn+ and Cr+, 
intermediate-ion transitions, A1++ and Fe ++ , and high -ion 
transitions, C 3+ and Si 3+ . To derive the intrinsic abundance 
measurements of the 11 detected elements, it is necessary to 
take into account both the ionization and dust depletion 
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Table 15. Spearman's tests and linear regressions 



Correlations for the whole sample of 
both the DLA and sub-DLA systems 


# of systems 


r 


a 




slope 


intercept 


log A(Si+) vs logA r (Fe+) 


39 


0.95 


3.2 x 10" 


-20 


1.04 ±0.04 


-0.10 ±0.58 


logiV(Al+) vs log A(Si+) 


19 


0.92 


3.4 x 10" 


-08 


0.97 ±0.08 


-0.97 ± 1.19 


logjV(Al+) vs logA(Fe+) 


17 


0.82 


4.9 x 10" 


-05 


0.89 ±0.13 


0.58 ± 1.83 


logiV(Al++)/A(Al+) vs logiV(H°) * 


30 


-0.67 


5.7 x 10- 


-05 


-0.59 ±0.13 


11.49 ± 2.84 



* Computed for the sample of DLAs only. 

References: This paper; Dessauges-Zavadsky et al. (2001a); Lopez et al. (1999,2002); Lu et al. (1995,1996b); 
Prochaska & Wolfe (1996,1997,1999); Prochaska et al. (2001); Peroux et al. (2002b) 



effects, since we are measuring gas-phase abundances. De- 
pletion is discussed in the wider context of Paper II. Here 
we analyse the ionization status of these systems. 

Several authors have investigated the ionization effects 
in the DLA systems (Viegas 1995; Howk & Sembach 1999; 
Izotov, Schaerer & Charbonnel 2001; Vladilo et al. 2001; 
Prochaska et al. 2002). Many of them were motivated by 
the observations of the intermediate-ion Al ++ , a tracer of 
moderately ionized gas, with the same velocity profiles as 
the low-ion transitions (Lu et al. 1996b; Prochaska & Wolfe 
1999; Prochaska et al. 2001), which suggest that ioniza- 
tion effects may be present in these systems. Indeed, to ex- 
plain the similarity of Al ill and other low-ionization species 
line profiles, Howk & Sembach (1999) and Izotov & Thuan 
(1999) proposed that these lines originate in the same ion- 
ized region or in a mix of neutral and ionized clouds, and 
stressed the importance of abundance corrections for ioniza- 
tion effects. Although the different approaches used to deal 
with the problem of photoionization in DLAs have led to 
slightly different conclusions, it is generally accepted that 
ionization corrections in DLAs are negligible, being below 
the measurement errors. 

In the present study, we are concerned about the ion- 
ization corrections for sub-DLA, i.e. systems with Hi col- 
umn densities lower than those of DLA systems, between 
10 19 < JV(Hl) < 2 x 10 20 cm" 2 . The check of ionization 
effects in these systems is very important, since their rela- 
tively low Hi column density implies that some of the gas 
might be ionized. 

4.1 A1 ++ /A1 + as Indicator of the Ionization State 

In six out of the 12 sub-DLA systems studied we detect the 
Aim A1854,1862 lines, and in four we also detect the Fein 
A1122 line. While neutral and singly ionized species with 
ionization potential < 13.6 eV can be present in the inner 
parts of H° regions, the production of Fe ++ and Al ++ re- 
quires photons with hv > 16.2 and 18.8 eV, respectively, 
which cannot so easily penetrate large Hi column densi- 
ties. The analysis of these intermediate-ion transition pro- 
files, however, shows that in all of the sub-DLAs but one 
at z a bs = 1.838 towards Q1101— 264 (see Section f3.2H . their 
velocity structure is almost indistinguishable from that of 
low-ion transition profiles, while the high-ion Civ and Siiv 
transitions have very different velocity profiles. This suggests 
that the Fe ++ and Al ++ and regions of neutral H in these 
systems are physically connected. 

The N{k\ ++ )/N(A\ + ) and N{¥e ++ ) / N(¥e + ) ratios 
can then be used as indicators of the ionization state. Vladilo 
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Figure 30. A1 ++ /A1 + column density ratios versus Hi column 
densities in sub-DLAs (filled symbols) and DLAs (open symbols) 
for all the systems found in the literature and studied in this 
paper. The squares correspond to the data points obtained with 
the A1+ column density measurements, while the circles and the 
triangles correspond to the data points computed from the Si + 
and Fe + column density measurements, respectively, according to 
the deduced linear regression parameters given in Table [T51 The 
dashed line represents a possible upper band limit of an increas- 
ing dispersion in the log Af(Al"'"+)/A r (Al+) ratios with lower Hi 
column densities. The slope of this upper band limit is equal to 
the slope of the observed log AT(Al+ + )/A r (Al+) versus log7V(H°) 
anti-correlation in DLAs. 



et al. (2001) studied the DLA ionization properties from 
the analysis of the A1 ++ /A1 + ratios in a data sample of 20 
DLAs, and found a log ; jV(Al ++ )/iV(Al + ) versus logAf(H°) 
anti-correlation. The existence of such an empirical anti- 
correlation implies that the neutral hydrogen column density 
is an indirect tracer of the ionization state of the gas. Con- 
sequently, if this anti-correlation is further observed in the 
sub-DLA Hi column density range, it would indicate that 
the sub-DLA systems require higher ionization corrections 
than absorbers in the DLA H I column density range. 

Therefore, we decided to investigate the extension of 
the different correlations found by Vladilo et al. (2001) in 
the DLA systems to the Hi column density range of the 
sub-DLA systems. For this purpose we have collected all 
the DLA and sub-DLA systems found in the literature for 
which measurements of the Al + and/or Al ++ column den- 
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sities exist and was performed on high resolution spectra 
obtained at the 8-10 m class telescopes. Most of them were 
derived from the HIRES/Keck spectra. We have then consid- 
ered only the systems which show similarities in their Al ++ 
and low-ionization species line profiles, since it is only in 
this case that the neutral and ionized clouds might be phys- 
ically related. Finally, we have adopted the values given by 
Vladilo et al. (2001) for the few cases of Al + column density 
measurements that they have revisited due to previously un- 
quoted saturations or contaminations by telluric lines. The 
results and the references for the adopted sample of both 
DLAs and sub-DLAs found in the literature and analysed 
in this paper are given in Table 1151 This table contains 
the Spearman's test correlation coefficients, r, the proba- 
bilities of null hypothesis (of no correlation), a, the slopes 
and the intercepts obtained for the linear regressions of 
log N(Si + ) versus log N(Fe + ), log N(A\ + ) versus log iV(Si + ), 
logAf(Al+) versus logA r (Fe+), and log iV(Al ++ )/iV(Al + ) 
versus logiV(H ). 

We confirm the existence of the logiV(Si + ) versus 
logAf(Fe+), logiV(Al + ) versus logiV(Si + ), and logiV(Al + ) 
versus log N(Fe + ) correlations found by Vladilo et al. (2001) 
in the DLA systems (see Table 1150 . These correlations are 
also satisfied by the sub-DLA systems. We then used the Si + 
(in priority, because the correlation with Al + has lower slope 
and intercept errors) or Fe + column densities as a proxy for 
iV(Al + ), when no iV(Al + ) measurement was available in the 
DLA or sub-DLA systems. In four sub-DLAs and several 
DLAs we took advantage of the correlations Si versus A1+ 
and Fe + versus Al + to indirectly estimate iV(Al + ) in order 
to increase the number of the log A r (Al ++ )/A r (Al + ) ratio 
measurements. 

Figure I3H1 shows the A1 ++ /A1 + column density ratios 
versus H° column densities both for the DLA (open sym- 
bols) and sub-DLA (filled symbols) systems. We confirm the 
log A r (Al ++ )/A r (Al+) versus logiV(H ) anti-correlation for 
the enlarged sample of DLA systems (30 objects instead of 
17 used by Vladilo et al. 2001) at the probability of null hy- 
pothesis of 5.7 x 10 -5 (correlation coefficient of r = —0.67; 
see Table ITSll . However, we show that at a lower Hi column 
density regime, this correlation is no longer valid. The sub- 
DLA systems have A1 ++ /A1 + ratios which are similar to 
those of the DLAs and not higher, as shown in Fig. 1301 with 
five measured values and seven upper limits (filled symbols). 
Only one system at z a b a = 3.565 towards PSS J2155+1358 
with logTV(Hl) = 19.37 ± 0.15 is really in agreement with 
the extrapolation of the DLA correlation to lower H I column 
densities (filled square in the upper left corner in Fig. 1301 . 

This implies that the neutral hydrogen column density 
is not a reliable indicator of the ionization state of the gas 
in absorption line systems as suggested by Vladilo et al. 
(2001), since it does not scale with the A1 ++ /A1 + column 
density ratio which is at a first approximation an indicator 
of the ionization level in the gas. On the contrary, Figure 1301 
shows that since the sub-DLA systems, despite their lower 
H I column densities, span approximately the same range of 
the A1 ++ /A1 + column density ratio values as the DLA sys- 
tems, they are likely to be affected by ionization corrections 
on average at the same level as the DLAs and not higher. 
However, we cannot claim that all of the absorption line sys- 
tems with Hi column densities between 10 19 and 2 x 10 20 
cm -2 (the sub-DLA H I column density definition) satisfy 



the same trend, as already shown by the example of the sub- 
DLA system at Zs , bs = 3.565 towards PSS J2155+1358. We 
can rather interpret the present data as showing an increas- 
ing dispersion in the log A'(Al ++ )/A r (Al + ) ratios with lower 
H° column densities. The plotted dashed line in Figure I3H1 
with the slope of the log iV(Al ++ )/jV(Ai + ) versus logiV(H ) 
anti-correlation observed in DLAs can be considered as an 
upper band limit of this increasing dispersion trend. If this 
trend is further confirmed by additional measurements, this 
would suggest that at high iV(H°) > 2 x 10 20 cm" 2 only low 
ionization effects are expected, while at lower iV(H°) be- 
tween 10 19 and 2 x 10 20 cm 2 one can find systems affected 
by low as well as high ionization corrections. 

4.2 Photoionization Models 

To investigate the ionization corrections in our sub-DLAs, 
we have used the CLOUDY software package (version 94.00, 
Ferland 1997), and we have computed photoionization mod- 
els assuming ionization equilibrium. In the models we have 
adopted a plane-parallel geometry for the gas cloud, the so- 
lar abundance pattern, and two different ionizing spectra, 
since the major uncertainty in determining the ionization 
effects in absorption line systems is the unknown shape of 
the ionizing spectrum. The two most likely origins for ion- 
izing UV photons are external background sources and in- 
ternal stellar sources — it is reasonable to expect some star 
formation in these systems. For the external ionizing source, 
we adopted the Haardt-Madau (hereafter HM) ionizing UV 
spectrum from background quasars with contributions from 
galaxies (Haardt & Madau 1996; Madau, Haardt & Rees 
1999). For the stellar ionizing source, we adopted the spec- 
trum of a local ionization which we consider to be a com- 
posite stellar spectrum plus a power law (hereafter SSP; see 
D'Odorico & Petitjean 2001). In the former case the ionizing 
photons have a harder spectrum. 

With each of these two ionizing spectra, we computed 
a series of CLOUDY models for each sub-DLA system in 
our sample by varying the H I column density of the plane- 
parallel 'cloud', the metallicity and the redshift according 
to the system considered, and by varying the volume den- 
sity of hydrogen, nn, through the variation of the ionization 
parameter, U, defined by 

TT _ $912 _ 471-J912 _ 5 J 9 i2/10~ 21 5 
U — — — — Z X 1U — , 

cuh ncnn nn /cm J 

where $912 is the surface flux of ionizing photons with hi; > 
1 Ryd, and J912 is the intensity of the incident radiation 
at 1 Ryd. In practice, we varied log U from —6 to 0, and 
obtained the theoretical column density predictions for all of 
the observed ions as a function of the ionization parameter. 

We then determined the ionization state of a given sub- 
DLA system by measuring adjacent ions of elements (e.g. 
A1++/A1+, Fe ++ /Fe+). Indeed, the value of the ionization 
parameter, U, at which the observed column density ratio 
of adjacent ions matches the predicted ratio, corresponds to 
the ionization parameter of the system studied. Therefore, 
the more numerous the measurements of column density ra- 
tios of adjacent ions of different elements, the more accurate 
is the determination of the ionization state of a system, all 
the more since the reliability of the derived photoioniza- 
tion models depends not only on the assumed hypothesis, 
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HM spectrum 
log N(HI) = 19.32; [Fe/H] = -0.84 



SSP spectrum 
log N(HI) = 19.32; [Fe/H] = -0.84 
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Figure 31. Photoionization models computed with the CLOUDY software package for the sub-DLA at z a b s = 2.557 towards Q1223+1753. 
This system is characterized by a relatively high metallicity [Fe/H] = — 0.84±0.15 and has one of the lowest H I column density (log N(H i) 
= 19.32 ±0.15) among the sub-DLA systems in our sample. On the left side, the models are the result obtained with the Haardt-Madau 
(HM) UV ionizing source, and on the right side, a stellar ionizing source (SSP) is assumed. In all panels the vertical lines indicate the 
observational constraint on U obtained from the A1++/A1+ (solid lines) and Fe++/Fe+ (long-dashed lines) ratios. The top panels show 
the predicted ionic column densities, 7V(X;), as a function of the ionization parameter, U. The horizontal dashed lines correspond to the 
observed column density values and the horizontal dotted lines correspond to the observed column density upper limits. The middle and 
bottom panels show the predicted abundance of low-ion Xi over Fe~*~ and H° relative to the intrinsic abundance of these two elements, 
as defined in Section 14.21 versus U, respectively. Departures of [Xj/Fe] and [X;/H] indicate that ionization corrections have to be added 
to/subtracted from the observed relative and absolute abundances, respectively. 



but also on the accuracy of the input atomic data for the 
CLOUDY code (e.g. dielectric recombination coefficients). 

Following the Prochaska et al. (2002) ionization analy- 
sis, we have further computed models giving the predicted 
abundance of low-ion Xi over H° and Fe + relative to the 
intrinsic abundance of these two elements, 



[X ; /H°] = log(JV(Xi)/JV(H°)) - log(X/H) , 
[Xi/Fe+] = log(JV(Xi)/JV(Fe + )) - log(X/Fe) , 

as a function of the ionization parameter, U. Departures of 
[Xi/H°] and [X;/Fe+] from zero indicate that ionization cor- 
rections have to be added to the observed absolute and rel- 
ative elemental abundances, respectively. When the correc- 
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Figure 32. Same as Fig. l31l for the sub-DLA at z a b s = 3.142 towards PSS J2155+1358. This system is characterized by a low metallicity 
[Fe/H] = —2.21 ± 0.21 and has one of the highest Hi column densities (log AT (Hi) = 19.94 ± 0.10) among the sub-DLA systems in our 
sample. In all panels the vertical line indicates the observational constraint on U obtained from the A1++/A1+ upper limit. 



tions to be applied are negative (positive), this implies that 
the measured [X/H] or [X/Fe] are underestimated (overes- 
timated) . The main advantage of these models compared to 
the models giving the predicted Xi column densities as a 
function of U is that they allow a direct quantification of 
the effects of photoionization on the absolute and relative 
abundances. 



4.3 Results: Ionization Corrections 

In the sample of the 12 sub-DLA systems and one bor- 
derline case between the DLA and sub-DLA systems stud- 
ied here, the Hi column densities range from log N (Hi) 



= 19.32 to 20.30 cm 2 , and the metallicities spread over 
1.5 dex (from [Fe/H] < -2.40 to [Fe/H] = -0.79). This large 
range of input parameters for the CLOUDY models allows 
us to highlight some general trends. Indeed, the computed 
CLOUDY models specific to each sub-DLA system in our 
sample show that for a given ionization parameter, U, the 
derived [Xi/H°] and [Xi/Fe + ] ionization corrections (as de- 
fined in Section r4.21 increase with lower H I column densities 
of systems for all low-ions Xi and for both the HM or SSP 
ionizing spectra, while no particular trend is observed with 
the metallicity of the systems. In addition, we notice that 
the ionization corrections are on the average more important 
in the models with a SSP radiation field than with a HM 
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spectrum for all of the systems studied. We show in Figs. 1311 
and !32l the CLOUDY models for two sub-DLA systems, one 
at z a bs = 2.557 towards Q1223+1753 with a very low H I 
column density logJV(Hl) = 19.32 ± 0.15 and a high [Fe/H] 
= —0.84 ± 0.15, and the second at z a b s = 3.142 towards 
PSS J2155+1358 with a high Hi column density logAT(Hi) 
= 19.94±0.10 and alow [Fe/H] = -2.21±0.21, respectively. 
They illustrate these general trends. 

Once we have computed the CLOUDY models, the de- 
termination of the ionization state is crucially dependent 
on the measurements of column density ratios of adjacent 
ions of elements, and thus on the detections of intermediate- 
ion transitions. In the sample of 13 systems studied, we 
were able to constrain the ionization parameter, U, for 
10 systems. However, we obtained for only four systems 
the iV(Al ++ )/iV(Ai + ) and/or AT(Fe ++ )/iV(Fe + ) measure- 
ments, while for the six remaining systems we derived only 
Ar(Al++)/iV(Al + ) and/or iV(Fe++)/iV(Fe + ) upper limits. 
The analysis of the ionization state of our sub-DLA sys- 
tems is thus severely limited by the fact that i) we do not 
have a measurement of a column density ratio of adjacent 
ions of an element for all of the sub-DLAs, and ii) for the 
majority of sub-DLAs we have only a measurement of the 
A1 ++ /A1 + ratio or even only an upper limit, and as has 
already been pointed out by Petitjean, Rauch & Carswell 
(1994), D'Odorico & Petitjean (2001) and Vladilo et al. 
(2001), the A1 ++ /A1 + ratio is not a very reliable indicator 
of the ionization parameter. Indeed, the recombination coef- 
ficient of Al + is likely overestimated (Nussbaumer & Storey 
1986), and in this case the predicted A1 ++ /A1 + ratio in the 
CLOUDY models is underestimated. This prevents us from 
obtaining a very accurate determination of the ionization 
parameter in our sub-DLA systems. In addition, with our 
set of data, we are not able to constrain the nature of the 
radiation field, since we cannot clearly differentiate which 
from the HM or SSP ionizing spectrum best represents the 
observations. 

In this context we choose to extract from our analysis 
general conclusions on the ionization corrections in the sub- 
DLA systems studied, rather than a discussion on individual 
systems which is limited by the partial information we have. 
Several important general remarks can be inferred from our 
study. First, as can be observed in Figs. I3ll and l32l the ion- 
ization corrections to apply to absolute and relative abun- 
dances are different ion to ion. Indeed, in the case of C + , O , 
S + or Cr + , for example, the ionization corrections [Xj/Fe + ] 
increase smoothly with the ionization parameter, U , and re- 
main relatively small for all of the U values, while for the 
Si + , Al + and Zn + ions they already become important at 
relatively low U values. Similarly, for the metallicities, the 
[Xj/H°] corrections are not significant for Fe + up to high U 
values, while they are much higher for Zn + at high U values. 
Secondly, for 9 systems out of 10 for which we were able to 
constrain the ionization parameter the ionization corrections 
are lower than 0.2 dex for all of the elements but Alii and 
Znll (in some cases; mainly if we assume that the SSP ra- 
diation field is the dominant ionizing source in our systems) 
down to Hi column densities of log N(H i) = 19.3 cm -2 (see 
Figs. EH and EU. In the case of Al + and Zn + the ionization 
corrections may be overestimated, since the Al + and Zn + 
recombination coefficients and ionization cross sections are 
rather uncertain (see Nussbaumer & Storey 1986 and Howk, 



Savage & Fabian 1999, respectively). Only in the sub-DLA 
system at z a b s = 3.565 towards PSS J2155+1358 are the 
ionization corrections important (filled square in the upper 
left corner in Fig. I3UII . For three sub-DLAs in our sample we 
could not constrain the ionization parameter, U. 



5 OVERVIEW OF THE CHEMICAL 
ABUNDANCES IN SUB-DLAS 

The constructed sample of 12 sub-DLA systems and one 
borderline case between the DLA and sub-DLA systems, 
and the homogeneous chemical abundance analysis of each 
system presented in this paper result in the first sub-DLA 
chemical abundance database ideal for the study of a num- 
ber of important issues. Indeed, 16 different ions of low-, 
intermediate- and high-ionization and 11 different elements 
- O, C, Si, N, S, Mg, Al, Fe, Ni, Zn and Cr - have been 
detected in the 13 systems studied, showing a variety of ele- 
ments as diverse as that observed in the DLAs. Tables [TBI a - ) 
and b) present a summary of the column density measure- 
ments of the low-ion transitions and the intermediate-/high- 
ion transitions, respectively, obtained in the 13 systems anal- 
ysed. These column densities were derived by summing the 
column densities of all of the detected velocity components 
in each ion as given in Tables I2T- I14I The cases in which we 
provide only lower limits correspond to the cases where the 
observed metal lines are saturated, and the cases in which 
we provide only upper limits correspond to cases where the 
observed metal lines are blended with Lya forest absorption 
lines. 4 a upper limits were derived in the cases where the 
metal lines are not detected and are identified by an asterisk. 

In Table lTTI we provide a summary of the absolute abun- 
dance measurements of the 12 sub-DLA systems and one 
borderline case between the DLA and sub-DLA systems, de- 
fined with respect to the solar meteoritic abundances from 
Grevesse & Sauval (1998). No ionization correction has been 
adopted since, as was shown in Section 01 for all of the 
systems studied for which we were able to constrain the 
ionization parameter (10 systems out of 13) except one at 
z a bs = 3.565 towards PSS J2155+1358, the ionization cor- 
rections are almost negligible, being at the level of the errors. 
Indeed, they are lower than 0.2 dex for all of the elements ex- 
cept Alii and Znll (in some cases, but likely overestimated 
because the Al + and Zn + recombination coefficients are very 
uncertain) down to Hi column densities of log iV(H i) = 19.3 
cm -2 . 

Fig. 1331 shows the distribution of iron and silicon abso- 
lute abundances, [Fe/H] and [Si/H], as a function of redshift 
and Hi column density for the 12 sub-DLAs in our sample. 
The studied sub-DLAs span a range of redshift from 1.8 to 
4.3 and of [Fe/H] metallicities from -0.79 to < -2.40, which 
is very similar to that observed in DLAs. We observe a trend 
of increasing metallicity with redshift, but no correlation be- 
tween metallicities and H I column densities is visible. An 
overall discussion on the chemical abundances of sub-DLAs, 
and a detailed comparison of the chemical, kinematic and 
clustering properties between the sub-DLAs and DLAs is 
given in Paper II. 

The lower H I column densities of the sub-DLA systems 
relative to the DLA systems give us the opportunity to have 
access to some elements which absorption lines are usually 
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Figure 33. The four plots show the distribution of the abun- 
dances of Fe and Si as a function of redshift and H I column den- 
sity, to illustrate the values sampled by the 12 sub-DLA systems 
in our sample. The overall abundances are presented in Table ITTI 



saturated in DLAs, but prevent us from measuring abun- 
dances of some other elements whose lines are already weak 
in DLAs. As examples of the former case, one can quote the 
carbon and oxygen. Indeed, C is almost not accessible in 
DLAs, since the Cn A1036,1334 lines are heavily saturated. 
Only very few tentative measurements have been obtained 
(Lopez et al. 1999,2002; Levshakov et al. 2002; Songaila & 
Cowie 2002). The O abundance is also very difficult to mea- 
sure in DLAs, since the most commonly detected O I line at 
A r cst = 1302 A is saturated and other O I lines are located 
in the Lya forest. On the contrary, in the sub-DLA systems 
one can more easily measure the abundances of these two 
important elements — we obtained four O and six C abun- 
dance measurements — which may provide insight into the 
understanding of the nature and the chemical evolution his- 
tory of these systems (see Paper II). On the other hand, the 
useful Zn abundance used as an indicator of the dust de- 
pletion in DLAs is very hard to measure in sub-DLAs given 
the weakness of the Znll lines. While numerous Zn abun- 
dance measurements exist in DLAs, we obtained only two 
Zn measurements in our sample of sub-DLAs. 



the Lya absorption line. Ten of them constitute a sample 
that we believe to be homogeneous and unbiased. 

• The detailed metal line analysis of each of the 12 sub- 
DLA systems and one borderline case between the DLA and 
sub-DLA systems resulted in the column density measure- 
ments of 16 ions of low-, intermediate- and high-ionization 
- Oi, Cn, Civ, Sin, Siiv, Ni, Sn, Mgi, Mgn, Alii, Aim, 
Fen, Fein, Nill, Znll and Cm. The detection of adjacent 
ions of the same element, Fe ++ /Fe + and A1 ++ /A1 + , allowed 
us to estimate the ionization parameters in 10 systems. Us- 
ing the CLOUDY software package, we computed photoion- 
ization models and found that for all of the 10 systems ex- 
cept one at z a b s = 3.565 towards PSS J2155+1358, the ion- 
ization corrections are almost negligible, being at the level 
of the errors. Indeed, they are lower than 0.2 dex for all of 
the elements except Alii and Znll (in some cases, but likely 
overestimated because the Al + and Zn + recombination co- 
efficients are very uncertain) down to Hi column densities 
of logTV(Hi) = 19.3 cm -2 . 

• We have obtained the first sub-DLA chemical abun- 
dance database which contains the abundance measure- 
ments of 11 different elements — O, C, Si, N, S, Mg, Al, 
Fe, Ni, Zn and Cr. The lower Hi column densities of sub- 
DLAs allowed us, in particular, to easily measure the O and 
C abundances, the O I and C II lines being usually saturated 
in DLAs. We obtained O abundance measurements in four 
systems and C in six systems. 

The constructed homogeneous and unbiased sample of 
the sub-DLA systems and the obtained sub-DLA chemical 
abundance database are ideal for the analysis of a number 
of important issues related both to the neutral gas mass and 
the metal content evolutions in the Universe. These issues 
are discussed in Paper II of this series, where we also present 
a detailed study of the sub-DLA kinematic, chemical and 
clustering properties, and their comparison with the prop- 
erties of the well studied DLA systems. We plan to extend 
the current sample of sub-DLA systems in order to increase 
the statistical weight of the obtained results both from the 
newly accessible data in the ESO UVES archives and from 
observations of pre-selected high-redshift quasars, since it is 
at high redshift z > 3.5 that the sub-DLAs are expected to 
contain a large fraction of the neutral hydrogen gas mass 
(Peroux et al. 2002a). 



6 CONCLUSION 

The findings by Peroux et al. (2003a) that the 'sub-DLA' 
systems with Hi column densities between 10 19 < A r (Hl) 
< 2 x 10 20 cm~ 2 contain a large fraction of the neutral 
hydrogen gas mass especially at z > 3.5, have motivated 
us to investigate systematically quasar spectra in the ESO 
UVES-VLT archives searching for these absorbers. We have 
obtained the following results: 

• Out of 22 high-resolution quasar spectra available to 
us in July 2001 in the ESO UVES-VLT archives and suit- 
able for the construction of a sample of sub-DLA systems, 
reduced and analysed with an identical approach, we have 
built the first sample of 12 sub-DLA systems in the redshift 
interval z = 1.8 — 4.3. These systems were identified in the 
quasar spectra both via an automated selection algorithm 
and visual criteria, like the presence of damping wings in 
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Table 16. a) Low- ion transition column density summary 



Name 


2 abs 


logJV(Hl) 


log JV(Ol) 


log JV(C II) 


logAT(Sill) 


log iV(Nl) 


logAT(Sn) 


logJV(Mgn) logiV(Aln) 


log iV(Fell) 


log iV(Ni II) 


log JV(Znll) 


log N(Cr II) 


tjllOl — 204 


1.838 


19.50(0.05) 


14.55(0.09) 


15.06(0.07) 


14.00(0.01) 




13.66(0. 11) 


14.07(0.03) 12.85(0.06) 


13.51 (0.02) 




< 11.39 


< 12.04 


^1223 + 1 f 53 


2.557 


19.32(0. 15) 


< 15.38 




14.38(0.02) 






13.20(0.03) 


13.98(0.03) 




< 11 .63 


< 11.96 


yi4uy+uyo 


2. (i()8 


19. 75(0. 10) 


IE it r t\ ii\ 
10.31(0.11) 


> 16.09 


1/1 1 Q / /i n/> \ 
14. 13(0 . DJ) 


< 13.49 


13.54(0.06) 




14.02(0. 13) 




< 11.34 


10 Iflln 1 E \ 

lz.4y(0. 15) 


1 A A A I m A 


2.087 


20. 18(0. 10) 


> 16.08 


> 16.30 


14.83(0.02) 






13.42(0.03) 


14.41(0.03) 








t^J14ol + lz3 




on 3rl'n i r \ 
ZU . 3U(0 . 1 5) 






> 15.46 








1 A 11 < r\ 

14.33(0.07) 


< 13.53 


11 Q E / n ii\ 
11.00(0.11) 


11 UD/n i/i\ 

12.00(0. 14) 


Q1451+123 


3.171 


19.70(0.15) 


< 14.80 


14.30(0.20) 


13.64(0.03) 






12.31(0.03) 


13.33(0.06) 




< 12.01* 


< 12.65* 


Q1511+090 


2.088 


19.47(0.10) 


< 15.30 




14.05(0.10) 






12.93(0.14) 


13.26(0.03) 








Q2059-360 


2.507 


20.21(0.10) 


15.54(0.21) 


14.87(0.19) 


13.73(0.11) 


12.66(0.18) 






13.47(0.03) 




< 11.54* 


< 12.19* 


Q2116-358 


1.996 


20.06(0.10) 


> 16.68 




15.33(0.12) 




< 14.95 


> 13.94 


14.77(0.09) 


13.57(0.06) 


12.29(0.09) 


12.92(0.13) 


PSS J2155+1358 


3.142 


19.94(0.10) 


< 14.71 


< 14.41 


13.65(0.07) 






12.33(0.10) 


13.23(0.19) 








PSS J2155+1358 


3.565 


19.37(0.15) 




14.34(0.12) 


13.66(0.04) 






12.16(0.12) 


< 12.47* 








PSS J2155+1358 


4.212 


19.61(0.10) 


14.50(0.05) 


13.95(0.06) 


13.25(0.04) 






11.92(0.17) 


12.93(0.23) 








PSS J2344+0342 6 ... 


3.882 


19.50(0.10) 




14.74(0.05) 


13.96(0.06) 






12.62(0.11) 


< 13.02* 
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tsa 

a. 

Co 



a The Mgl column density is also measured in this sub-DLA: log iV(Mgl) 
6 If the components 1—6 do not belong to the sub-DLA (see Section 13, 111 . 



= 11.89 ± 0.02. 
then: logJV(Cn) 



14.55 ± 0.04, log N(Si n) = 13.70 ± 0.04 and log iV(Al n) = 12.36 ± 0.07. 



b) Intermediate- and high-ion transition column density summary 



Name 


z abs 


log JV(Hl) 


logJV(Alm) 


log JV(Fe m) 


log JV(C IV) 


logiV(Silv) 


Q1101-264 


1.838 


19.50(0.05) 


12.34(0.04) 




14.22(0.04) 


13.84(0.05) 


Q1223+1753 


2.557 


19.32(0.15) 


12.32(0.13) 


14.31(0.13) 


14.44(0.04) 


13.99(0.02) 


Q1409+095 


2.668 


19.75(0.10) 


11.89(0.05) 


13.17(0.12) 






Q1444+014 


2.087 


20.18(0.10) 


< 12.06* 


< 13.90* 


13.22(0.05) 


12.94(0.04) 


Q1451+123 


2.255 


20.30(0.15) 


12.70(0.05) 








Q1451+123 


3.171 


19.70(0.15) 


< 11.57* 




13.48(0.13) 




Q1511+090 


2.088 


19.47(0.10) 


< 11.84* 


< 13.75* 






Q2059-360 


2.507 


20.21(0.10) 


< 11.89* 








Q2116-358 


1.996 


20.06(0.10) 


13.25(0.09) 


< 14.60 


> 15.28 


14.38(0.02) 


PSS J2155+1358 


3.142 


19.94(0.10) 


< 11.57 








PSS J2155+1358 


3.565 


19.37(0.15) 


12.71(0.09) 


< 13.81 


14.38(0.10) 


14.05(0.04) 


PSS J2155+1358 


4.212 


19.61(0.10) 








12.93(0.12) 


PSS J2344+0342 


3.882 


19.50(0.10) 


< 12.32* 




14.19(0.05) 


13.64(0.07) 



4 a upper limit (corresponding to non-detections). 



Table 17. Absolute abundance summary 



Name 


z abs 


log JV(Hl) 


[O/H] 


[C/H] 


[Si/H] 


[N/H] 


[S/H] 


[Mg/H] 


[Al/H] 


[Fc/H] 


[Ni/H] [Zn/H] 


[Gr/H] 


Q1101-264 


1.838 


19.50(0.05) 


-1.78(0.12) 


-0.96(0.10) 


-1.06(0.05) 




-1.04(0.13) 


-1.01(0.06) 


-1.14(0.08) 


-1.49(0.05) 


< -0.78* 


< -1.15* 


Q1223+1753 


2.557 


19.32(0.15) 


< -0.77 




-0.50(0.15) 








-0.61(0.15) 


-0.84(0.15) 


< -0.36* 


< -1.05* 


Q1409+095 


2.668 


19.75(0.10) 


-1.27(0.16) 


> -0.18 


-1.18(0.10) 


< -2.18 


-1.41(0.13) 






-1.23(0.16) 


< -1.08* 


-0.95(0.18) 


Q1444+014 


2.087 


20.18(0.10) 


> -0.93 


> -0.40 


-0.91(0.10) 








-1.25(0.10) 


-1.27(0.10) 






Q1451 + 123 


2.255 


20.30(0.15) 






> -0.40 










-1.47(0.17) 


< -1.02 -1.12(0.19) 


-1.11(0.20) 


Q1451+123 


3.171 


19.70(0.15) 


< -1.73 


-1.92(0.26) 


-1.62(0.15) 








-1.88(0.15) 


-1.87(0.16) 


< -0.36* 


< -0.74* 


Q1511+090 


2.088 


19.47(0.10) 


< -1.00 




-0.98(0.14) 








-1.03(0.17) 


-1.71(0.10) 






Q2059-360 


2.507 


20.21(0.10) 


-1.50(0.24) 


-1.86(0.22) 


-2.04(0.15) 


-3.47(0.21) 








-2.24(0.10) 


< -1.34* 


< -1.71* 


Q2116-358 


1.996 


20.06(0.10) 


> -0.21 




-0.29(0.16) 




< -0.31 




> -0.61 


-0.79(0.13) 


—0.74(0.12) —0.44(0.14) 


-0.83(0.16) 


PSS J2155+1358 


3.142 


19.94(0.10) 


< -2.06 


< -2.05 


-1.85(0.12) 








-2.10(0.14) 


-2.21(0.21) 






PSS J2155+1358 c ... 


3.565 


19.37(0.15) 




— 1.55(0.20) 


-1.27(0.16) 








-1.70(0.19) 


< -2.40* 






PSS J2155+1358 


4.212 


19.61(0.10) 


-1.94(0.13) 


-2.18(0.13) 


-1.92(0.11) 








-2.18(0.20) 


-2.18(0.25) 






PSS J2344+0342 


3.882 


19.50(0.10) 




-1.28(0.13) 


-1.10(0.12) 








-1.37(0.15) 


< -1.98* 







* 4 a upper limit (corresponding to non-detections). 

c This system is very likely affected by high ionization corrections (see Section HI. 
Abundances relative to the solar values of Grevesse & Sauval (1998). 
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